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1.  SUMMARY 

This  is  the  first  Annual  Report  under  ONR  Contract  N00014-84-0468,  work  unit  number 
NR064-725  conducted  by  Drs.  H.  Murakami  and  G.  A.  Hegemier  at  the  University  of  Califor¬ 
nia.  San  Diego,  during  the  period  from  July  1,  1984  to  June  30,  1985.  The  cognizant  ONR  pro¬ 
gram  manager  is  Dr.  A.  S.  Kushner. 

he  research  is  conducted  in  response  to  the  need  for  microstructural-based  theories  which 
furnish  increased  simulation  capabilities  of  metal-matrix  composite  structures  in  both  static  and 
dynamic  regimes  with  a  minimum  of  model  parameters  to  be  experimentally  determined.  The 
objective  of  the  research  is  to  develop  a  nonlinear  model  of  binary  metal-matrix  composites  that 
will  provide  greater  accuracy  than  existing  models  in  linear  and  nonlinear  regimes  for  static  and 
dynamic  loading.  Furthermore,  in  an  effort  to  accommodate  practical  structural  configurations 
a  plate  theory  is  developed  for  fiber-reinforced  composite  plates. 

The  progress  made  during  the  first  year’s  effort  toward  achieving  the  above  objectives 
includes:  |fl)  the  development  and  validation  of  a  linear  mixture  model  which  accounts  for  effec¬ 
tive  moduli  and  harmonic  wave  dispersion,  |2)  testing  of  methodologies  for  including  nonlinear 
material  responses,  and  ^3)  a  significant  improvement  of  existing  laminated  composite  theories 
to  account  for  inelastic  responses  driven  by  in-plane  strains. 


'T 


2.  RESEARCH  OBJECTIVE 

The  ultimate  objective  of  the  research  described  here  is  to  construct  an  advanced,  non¬ 
linear  nonphenomenological  model  of  binary  metal-matrix  composites  that  will  provide  greater 
accuracy  than  existing  models  in  linear  and  nonlinear  response  regimes  for  static  and  dynamic 
loading.  The  term  "nonphenomenological’  implies  a  model  that  is  capable  of  synthesizing  the 
global  properties  of  composites  from  a  knowledge  of  the  matrix  and  fiber  properties,  the  fiber- 
matrix  interface  properties,  and  the  geometry  of  the  fiber  reinforcement.  In  addition  to  the 
development  of  a  continuum  model,  a  plate  theory  is  developed  for  laminated  composites  in 
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which  each  layer  consists  of  a  unidirectionally  fiber-reinforced  metal-matrix  composite  laminae. 
This  effort  is  to  improve  the  simulation  capabilities  of  metal-matrix  composite  structures  by 
introducing  the  advanced  constitutive  model  to  laminated  composite  plates. 

The  specific  research  objectives  of  the  work  performed  during  the  present  reporting  period 
were  as  follows: 

1.  Construction  of  a  linear  mixture  model  -  Develop  a  dispersive  linear  mixture  model 
for  arbitrary  wave  motion,  and  perform  a  validation  study  of  the  model. 

2.  Development  of  a  nonlinear  mixture  model  -  Explore  an  efficient  methodology  for 
including  inelastic  responses  of  the  composite,  such  as  plastic  deformation  of  the 
matrix,  debonding  and  slip  at  the  fiber-matrix  interface,  and  fiber  breakage. 

3.  Improvement  of  existing  laminated  composite  plate  theories  -  Develop  an  improved 
laminated  composite  plate  theory  which  can  accommodate  the  above  inelastic  consti¬ 
tutive  model  for  metal-matrix  composite  laminae. 


S.  CURRENT  STATUS  OF  RESEARCH 

The  progress  made  during  the  report  period  toward  achieving  the  research  objectives 
described  earlier  is  summarized  in  this  section.  First,  the  basic  technical  approach  being  fol¬ 
lowed  to  meet  the  objectives  is  outlined.  Then,  the  progress  made  toward  developing  and  vali¬ 
dating  an  advanced  mixture  theory  with  microstructure  for  metal-matrix  composites  is 
described.  Next,  the  effort  to  develop  a  laminated  composite  plate  theory  with  improved  in¬ 
plane  responses  for  use  in  connection  with  the  above  constitutive  model  is  summarized. 

S.l  Approach 

The  nonlinear  response  of  metal-matrix  composites  is  largely  dominated  by  complex 
interaction  between  the  fiber  and  the  metal-matrix.  Consequently,  an  accurate  model  of  metal- 
matrix  composites  must  be  capable  of  accounting  for  such  interactions.  Further,  in  an  effort  to 
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minimize  the  number  and  types  of  tests  necessary  to  define  the  parameters  of  a  given  model,  it 
is  highly  desirable  that  it  be  nonphenomenological,  i.e.,  that  the  global  properties  of  the  compo¬ 
site  be  synthesized  from  the  constitutive  properties  of  the  fiber  and  matrix,  the  fiber-matrix 
interface  conditions,  and  the  geometry  of  fiber  reinforcement.  A  candidate  modeling  approach 
that  satisfies  the  above  objective  is  the  "mixture  theory  with  microstructure".  According  to  the 
mixture  concept,  the  fiber  and  matrix  are  modeled  at  each  instant  of  time  as  superposed  con- 
tinua  in  space.  Each  continuum  is  allowed  to  undergo  individual  deformations.  The  micros¬ 
tructure  of  an  actual  composite  is  then  simulated  by  specifying  the  nature  of  the  interaction 
between  the  continua.  The  key  element  for  the  development  of  a  mixture  model  for  composites 
with  periodic  microstructure  is  an  asymptotic  procedure  called  'Vnultivariable  asymptotic  expan¬ 
sions".  This  mathematical  technique  with  a  'Smoothing"  operation,  leads  to  the  desired  mixture 
forms. 

An  improved  laminated  composite  plate  theory  that  can  simulate  in-plane  responses  accu¬ 
rately  is  developed  by  introducing  a  new  displacement  microstructure  over  the  thickness  of 
plates  and  by  using  Reissner’s  new  mixed  variational  principle  (Reissner,  1984)  which  automati¬ 
cally  yields  the  shear  correction  factors  of  shear  deformable  plate  theories  (Reissner,  1946,  Mind- 
lin,  1951). 


S.2  Development  of  Linear  and  Nonlinear  Mixture  Models  for  Metal-Matrix  Composites 

For  fibrous  composites,  wave  dispersion  has  been  amply  demonstrated  via  ultrasonic  tech¬ 
niques  by  such  investigators  as  Tauchert  and  Guzelse  (1972),  and  Sutherland  and  Lingle  (1972). 

Simulation  of  response  phenomena  associated  with  the  material  microstructure,  such  as 
wave  dispersion,  requires  a  higher-order  continuum  description.  Several  such  models  have  been 
proposed,  some  phenomenological,  some  nonphenomenological. 

A  higher-order  continuum  model  which  simulates  wave  dispersion  was  first  proposed  by 
Achenback  and  Herrmann  (1968)  for  unidirectionally  fiber-reinforced  composites.  This  theory, 
called  the  'Effective  stiffness  theory",  has  been  further  studied  and  applied  to  fibrous  composites 


by  Bartholomew  and  Torvick  (1972),  Hlavacek  (1975),  Achenback  (1976),  and  Aboudi  (1981). 
The  aforementioned  work  concerned  linear  materials.  By  modifying  the  original  methodology, 
Aboudi  (1982,  1985)  extended  the  linear  model  to  account  for  inelastic  responses  of  the  compo¬ 
site  constituents. 

In  addition  to  the  effective  stiffness  modeling  concept,  a  mixture  approach  has  been  fol¬ 
lowed  by  a  number  of  investigators.  A  phenomenological  version  of  this  model  type  was 
adopted  by  Martin,  Bedford  and  Stern  (1971).  Deterministic,  nonphenomenologica)  mixture 
theories  were  introduced  by  Hegemier,  Gurtman  and  Nayfeh  (1973),  Hegemier  and  Gurtman 
(1974),  and  Murakami,  Maewal  and  Hegemier  (1979).  Although  capable  of  simulating  nonlinear 
component  responses  and  interfacial  slip,  this  work  was  limited  to  waveguide-type  problems. 
This  limitation  was  removed  in  the  mixture  theory  developed  for  laminated  composites  by 
Hegemier,  Murakami  and  Maewal  (1979),  and  Murakami,  Maewal  and  Hegemier  (1982).  In 
their  papers,  it  was  demonstrated  that  the  mixture-type  model  was  capable  of  simulating  har¬ 
monic  wave  dispersion  in  laminated  composites  more  accurately  than  the  effective  stiffness 
theories.  Further,  the  mixture-type  model  requires  fewer  governing  equations.  The  accuracy 
and  efficiency  of  the  mixture  theory  is  due  to  the  use  of  appropriate  displacement  and  stress 
microstructural  fields,  and  a  judicious  smoothing  technique.  These  are  obtained  by  an  asymp¬ 
totic  procedure  with  multiple  scales.  This  procedure  yields  a  series  of  microboundary  value 
problems  (MBVP’s)  defined  over  a  unit  cell,  which  in  term  represents  the  (periodic)  microstruc¬ 
ture  of  a  composite.  The  lowest  order  version  of  the  MBVP  method  is  equivalent  to  the  '0(1) 
homogenization  theory”  summarized  by  Bensoussan,  Lions,  and  Papanicolaou  (1978),  and 
Sanchez-Palencia  (1980).  The  latter,  while  it  generates  appropriate  static  moduli,  is  nondisper- 
sive.  Simulation  of  wave  dispersion  requires  at  least  a  theory  which  is  classified  as  an  0(e) 
homogenization  theory  in  which  e  denotes  the  representative  ratio  of  micro-to-microdimensions 
of  a  composite. 

To  date  an  O(e)  mixture  theory  has  not  been  constructed  for  fibrous  composites  subject  to 
arbitrary  wave  motion.  Construction  and  validation  of  such  a  3D  model  for  unidirectional 
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binary  composites  with  periodic  microstructure  is  the  objective  of  the  research.  To  facilitate 
this  task,  the  asymptotic  procedure  with  multiple  scales  noted  previously  is  combined  with  a 
variational  technique  (Murakami,  1985).  Following  development  of  the  basic  equations,  the 
dispersion  of  time-harmonic  waves  is  studied  and  the  results  are  compared  with  experimental 
data  for  boron/epoxy  (Tauchert  and  Guzelse,  1972)  and  tungusten/aluminum  (Sutherland  and 
Lingle,  1972)  composites.  The  good  correlation  obtained  with  experimental  data  indicates  that 
the  proposed  mixture  model  furnishes  a  basic  tool  by  which  dynamic  responses  of  elastic  compo¬ 
sites  can  be  investigated. 

Following  the  development  of  a  linear  mixture  model  for  metal-matrix  composites,  an 
extension  of  the  model  to  include  material  nonlinearities  has  been  attempted.  Future  work  will 
include  completion  of  development  of  a  nonlinear  mixture  model  with  validation  studies.  As 
part  of  this  effort,  the  constraint  hardening  and  fiber  breakage  will  be  incorporated. 

S.S  Development  of  An  Improved  Laminated  Composite  Plate  Theory 

The  application  of  metal-matrix  composites  in  the  form  of  laminated  plates  has  created  a 
demand  for  the  development  of  a  laminated  composite  plate  theory  in  which  each  layer  may 
experience  plastic  deformation  with  constraint  hardening  (Dvorak  and  coworkers,  1976,  1984), 
and  transverse  cracking.  In  order  to  simulate  the  inelastic  response  of  each  layer,  plate  theories 
should  be  capable  of  predicting  accurately  in-plane  strains  which  yield  inelastic  responses. 

In  a  series  of  papers,  Pagano  ( 1970a, b)  derived  exact  elasticity  solutions  for  bidirectional 
composites  for  the  problems  of  cylindrical  bending  and  simply  supported  rectangular  plates. 
Pagano  showed  the  importance  of  the  transverse  shear  effect  for  the  predictions  of  accurate 
plate  deflections  and  the  necessity  of  improving  assumptions  for  in-plane  displacements,  which 
are  assumed  to  be  linear  across  the  thickness  of  the  plate  in  the  Kirchhoff  as  well  as  the 
Reissner-Mindlin  shear  deformable  plate  theories.  Since  the  development  of  laminated  plate 
theories,  including  the  effect  of  the  transverse  shear  by  Yang,  Norris  and  Stavsky  (1966)  and 
Whitney  and  Pagano  (1970),  many  higher  order  laminated  plate  theories  have  been  proposed. 
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Historical  accounts  of  such  efforts  may  be  found  in  the  articles  by  Seide  (1980),  Bert  (1984),  and 
Reddy  (1984).  However,  only  a  few  attempts  have  been  made  to  improve  the  in-plane  strain 
responses. 

New  high-order  plate  theories  have  been  developed  with  the  help  of  a  new  variational  prin¬ 
ciple  (Reissner,  1984).  The  improvement  of  the  in-plane  responses  is  achieved  by  including  a 
zigzag  shaped  C°  function  to  approximate  the  thickness  variation  of  the  in-plane  displacements. 

In  any  approximate  plate  theory,  kinematic  assumptions  require  corresponding  constitutive 
assumptions  for  transverse  stresses.  As  an  example,  both  the  KirchhofT  and  the  Reissner- 
Mindlin  plate  theories  adopt  a  displacement  field  which  satisfies  a  state  of  plane  strain  in  the 
thickness  direction.  In  the  Kirchhoff  theory  the  stress  field  is  assumed  to  be  in  a  state  of  plane 
stress  which  means  that  all  transverse  stresses  are  zero.  In  the  Reissner-Mindlin  plate  theory 
only  the  transverse  normal  stress  is  set  to  be  zero  in  the  stress-strain  relations.  This,  in  turn, 
implies  that  Hook’s  law,  as  it  is,  cannot  be  used.  As  a  result,  a  more  complicated  displacement 
assumption  over  the  thickness  of  the  plate  must  be  introduced  with  suitable  stress  assumptions 
which  are  not  trivial.  Reissner’s  new  mixed  variational  principle  (1984)  provides  a  solution  to 
the  above  shortcoming;  it  is  a  variational  principle  for  arbitrary  displacement  and  transverse 
stresses,  in  which  the  original  three  dimensional  stress-strain  relations  can  be  used. 

By  using  the  two  key  elements:  a  zigzag  shaped  C°  interpolation  function  for  the  displace¬ 
ment  variation  over  the  thickness  of  the  plate  and  the  application  of  Reissner’s  new  mixed  vari¬ 
ational  principle,  new  shear  deformable  laminated  plate  theories  with  improved  in-plane 
responses  were  developed  (Murakami,  1985,  Toledano  and  Murakami,  1985).  Copies  of  these 
references  are  included  in  Appendix  B,  together  with  the  paper  which  compared  the  difference  of 
the  two  plate  theories  of  different  order  (Murakami  and  Toledano,  1985). 


The  overall  thrust  and  implication  of  the  above  work  is  that  by  the  new  laminated  compo¬ 
site  plate  theories  it  is  now  possible  to  carry  out  nonlinear  plate  analyses  in  which  some  layers 
experience  plastic  deformations.  For  the  composite  plates  made  of  metal-matrix  composite  lam¬ 
inae  the  new  plate  theories  can  accommodate  the  new  constitutive  model  developed  by  the  non- 
phenomenological  mixture  theory. 


4.  PUBLICATIONS 

The  following  papers  were  prepared,  and  submitted  for  publication,  during  the  reporting 
period  covered  by  this  report: 

1.  Murakami,  H.  and  G.  A.  Hegemier,  "A  Mixture  Model  for  Unidirectionally  Fiber- 
Reinforced  Composites,”  submitted  for  publication. 

2.  Murakami,  H.,  "Laminated  Composite  Plate  Theory  with  Improved  In-Plane 
Responses,”  Proceedings  of  the  1985  PVP  Conference,  ASME,  PVP  Vol.  98-2,  1985, 
pp.  257-263,  also  submitted  to  ASME  Journal  of  Applied  Mechanics. 

3.  Toledano,  A.  and  H.  Murakami,  ”A  High-Order  Laminated  Plate  Theory  with 
Improved  In-Plane  Responses,”  submitted  for  publication. 

The  following  paper  was  prepared  for  a  cancelled  symposium  during  the  reporting 
period  and  will  be  submitted  to  an  appropriate  symposium. 

4.  Murakami,  H.  and  Toledano.  A.,  ”An  Improved  Laminated  Composite  Plate  Theory," 


will  be  submitted. 


5.  INTERACTIONS 


The  following  is  a  list  of  the  presentations  at  meetings  and  conferences  by  the  principal 
investigators  which  occurred  during  the  reporting  period  on  issues  related  to  the  research  done 
under  the  present  contract: 

1.  Murakami,  H.,  "A  Mixture  Model  for  Metal-Matrix  Composites,"  oral  presentation  at 
the  Tenth  Annual  Mechanics  of  Composite  Review,  Dayton,  Ohio,  October  15-17, 
1984. 

2.  Murakami,  H.,  'Some  Basic  Inelastic  Response  Features  of  the  New  Endochronic 
Theory,"  oral  presentation  at  the  21st  Annual  Meeting  of  the  Societv  of  Engineering 
Science,  VPI,  Blacksburg,  Virginia,  October  15,  1984. 

3.  Murakami,  H.,  "Laminated  Composite  Plate  Theory  with  Improved  In-Plane 
Responses,"  oral  presentation  at  the  1985  PVP  conference,  ASME,  New  Orleans, 
Louisiana,  June  24-26,  1985. 

6.  LIST  OF  PROFESSIONAL  PERSONNEL 

Scientific  personnel  supported  by  the  contract  during  the  reporting  period  are 

1.  !  .incipal  Investigators:  Dr.  H.  Murakami,  Assistant  Professor  of  Applied  Mechanics, 
and  Dr.  G.  A.  Hegemier,  Professor  of  Applied  Mechanics. 

2.  Research  Assistants:  Mr.  Akira  Akiyama,  Mr.  Albert  Toledano,  and  Mr.  Thomas 
Impelso,  PhD  students  in  Engineering  Sciences  (Applied  Mechanics). 
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ABSTRACT 


A  binary  mixture  theory  with  microstructure  is  constructed  Tor  unidirectionally  fiber-reinforced 
elastic  composites.  Model  construction  is  based  on  an  asymptotic  scheme  with  multiple  scales  and  the 
application  of  Reissner's  new  mixed  variational  principle  (1984).  In  order  to  assess  the  accuracy  of  the 
model,  comparison  of  the  mixture  model  predictions  with  available  experimental  data  on  dispersion  of 
harmonic  waves  is  made  for  boron/epoxy  and  tungusten/aluminum  composites.  Formulas  for  the 
effective  moduli  are  also  presented,  and  the  results  are  compared  with  test  data  and  other  available 
predictions. 
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1.  Introduction 

With  the  advent  of  high  strength  and  stiffness  fibers  such  as  boron  and  carbon,  and  the  develop¬ 
ment  of  techniques  for  binding  such  materials  to  plastic  or  metal,  fibrous  composites  have  become 
important  elements  of  modern  structures.  Such  composites,  due  to  their  microstructural  heterogeneity, 
may  exhibit  response  phenomena  for  some  environments  that  are  not  observed  for  homogeneous 
materials.  An  example  of  these  phenomena  for  dynamic  environments  is  wave  dispersion,  and  under¬ 
standing  of  which  is  important  both  from  the  standpoints  of  direct  response  prediction  and  indirect  ana¬ 
lyses  associated  with  such  topics  as  nondestructive  testing.  For  fibrous  composites,  wave  dispersion  has 
been  amply  demonstrated  via  ultrasonic  techniques  by  such  investigators  as  Tauchert  and  Guzelse 
(1972),  and  Sutherland  and  Lingle  (1972). 

Simulation  of  response  phenomena  associated  with  the  material  microstructure,  such  as  wave 
dispersion,  requires  a  higher-order  continuum  description.  Several  such  models  have  been  proposed, 
some  phenomenological,  some  nonphenomenological. 

A  higher-order  continuum  model  which  simulates  wave  dispersion  was  first  proposed  by  Achen- 
back  and  Herrmann  (1968)  for  unidirectionally  fiber-reinforced  composites.  This  theory,  called  the 
"effective  stiffness  theory,  has  been  further  studied  and  applied  to  fibrous  composites  by  Bartholomew 
and  Torvick  (1972),  Hlavacek  (1975),  Achenback  (1976),  and  Aboudi  (1981).  The  aforementioned 
work  concerned  linear  materials.  By  modifying  the  original  methodology,  Aboudi  (1982,  1983) 
extended  the  linear  model  to  account  for  inelastic  responses  of  the  composite  constituents. 

In  addition  to  the  effective  stiffness  modeling  concept,  a  mixture  approach  has  been  followed  by  a 
number  of  investigators.  A  phenomenological  version  of  this  model  type  was  adopted  by  Martin,  Bed¬ 
ford  and  Stern  (1971).  Deterministic,  nonphenomenological  mixture  theories  were  introduced  by 
Hegemier,  Gurtman  and  Nayfeh  (1973),  Hegemier  and  Gurtman  (1974),  ar.  Murakami,  Maewal  and 
Hegemier  (1979).  Although  capable  of  simulating  nonlinear  component  responses  and  interfacial  slip, 
this  work  was  limited  to  waveguide-type  problems.  This  limitation  was  removed  in  the  mixture  theory 
developed  for  laminated  composites  by  Hegemier,  Murakami  and  Maewal  (1979),  and  Murakami, 
Maewal  and  Hegemier  (1982).  In  their  papers,  it  was  demonstrated  that  the  mixture-type  model  was 
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capable  of  simulating  harmonic  wave  dispersion  in  laminated  composites  more  accurately  than  the 
effective  stiffness  theories.  Further,  the  mixture-type  model  requires  fewer  governing  equations.  The 
accuracy  and  efficiency  of  the  mixture  theory  is  due  to  the  use  of  appropriate  displacement  and  stress 
microstructural  fields,  and  a  judicious  smoothing  technique.  These  are  obtained  by  an  asymptotic  pro¬ 
cedure  with  multiple  scales.  This  procedure  yields  a  series  of  microboundary  value  problems  (MBVP’s) 
defined  over  a  unit  cell,  which  in  turn  represents  the  (periodic)  microstructure  of  a  composite.  The 
lowest  order  version  of  the  MBVP  method  is  equivalent  to  the  "0(1)  homogenization  theory  summar¬ 
ized  by  Bensoussan,  Lions,  and  Papanicolaou  (1978),  and  Sanchez-Palencia  (1980).  The  latter,  while  it 
generates  appropriate  static  moduli,  is  nondispersive.  Simulation  of  wave  dispersion  requires  at  least  a 
theory  which  is  classified  as  an  0(e)  homogenization  theory  in  which  e  denotes  the  representative  ratio 
of  micro-to-macrodimensions  of  a  composite. 

To  date  an  0(e)  mixture  theory  has  not  been  constructed  for  fibrous  composites  subject  to  arbi¬ 
trary  wave  motion.  Construction  and  validation  of  such  a  3D  model  for  unidirectional  binary  compo¬ 
sites  with  periodic  microstructure  are  the  objective  of  this  paper.  To  facilitate  this  task,  the  asymptotic 
procedure  with  multiple  scales  noted  previously  is  combined,  with  a  variational  technique  (Murakami, 
1985).  Following  development  of  the  basic  equations,  the  dispersion  of  time-harmonic  waves  is  studied 
and  the  results  are  compared  with  experimental  data  for  boron/epoxy  (Tauchert  and  Guzelse,  1972) 
and  tungusten/aluminum  (Sutherland  and  Lingle,  1972)  composites.  The  good  correlation  obtained 
with  experimental  data  indicates  that  the  proposed  mixture  model  furnishes  a  basic  tool  by  which 
dynamic  responses  of  elastic  composites  can  be  investigated.  While  the  model  construction  procedure 
is  applicable  to  inelastic  component  response  and  interface  slip,  extension  and  investigation  of  the  non¬ 
linear  problem  is  deferred  to  later  publications. 

2.  Formulation 

Consider  a  domain  V  which  contains  a  uniaxial  periodic  array  of  fibers  embedded  in  the  matrix,  as 
shown  in  Fig.  1.  Let  a  rectangular  reference  system  3ft,  x 2,  x3  be  selected  with  jc,  in  the  axial  direction 
of  the  fibers.  In  the  x2,  x3 -plane,  a  typical  cell  that  represents  the  geometrical  microstructure  of  the 


composite  is  shown  in  Fig.  2  for  a  hexogonal  array. 


For  notational  convenience  forms  (  )<a),  a  *  1,2  denote  quantities  associated  with  material  a 
with  a  -  1  representing  fiber  and  a  -  2  matrix.  Cartesian  indicia!  notation  will  be  employed  in  which 
Latin  indices  range  from  1  to  3  and  repeated  indices  imply  the  summation  convention  unless  otherwise 
stated.  In  addition,  the  notations  C),  =  8  C  )/8x,  and  C)j  =  80/87  will  be  employed  in  which  7 
represents  time.  Quantities  of  the  form  (”)  and  (  )  denote  dimensional  and  nondimensional  variables, 
respectively. 

The  governing  relations  for  the  displacement  vector  «,<o)  and  the  stress  tensor  cr^’  in  the  two 
constituents  are: 


(a)  Equations  of  motion 


j  P  uui 


where  is  the  mass  density; 

(b)  Constitutive  relations 


(1) 


e£>  +  ef'  (2) 

wheu  \(a),  fl(a)  are  Lame's  constants,  e,Ja>  is  the  infinitesimal  Cauchy  strain,  and  80  is  the  Kronecker 
delta; 

(c)  Strain-displacement  relations 

\  my  +  Tiff)  ;  (3) 


(d)  Interface  continuity  relations 

ty(,)  -  m,<2)  ,  ^(,V»  -  on  3 

where  vj11  =  0  on  the  fiber-matrix  interface  3  '•> 

(e)  Initial  conditions  at  7—  0  and  appropriate  boundary  data  along  the  boundary  8  V . 


(4) 


Conditions  (a)  •  (e)  define  a  well  posed  initial  boundary  value  problem.  However,  due  to  the 
large  number  of  fiber-matrix  interfaces  the  direct  solution  to  this  problem  is  extremely  difficult.  The 
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objective  of  the  subsequent  analysis  is  to  alleviate  such  difficulties  by  deriving  a  set  of  partial  differential 


equations  with  constant  coefficients  whose  solution  can  be  utilized  to  approximate  the  solution  of  the 


problem.  To  this  end,  it  will  be  convenient  to  nondimensionalize  the  basic  equations  by  using  the  fol¬ 


lowing  quantities: 


typical  macrosignal  wavelength 


typical  fiber  spacing  or  cell  dimension 


C(m)  ,  P(m) 


reference  wave  velocity  and  macrodensity 


E(m)  =  P(m)  C}m)  reference  modulus 


Um)  =  A JCon) 


typical  macrosignal  travel  time 


€  =  A/A 


ratio  of  micro-to-macrodimensions. 


With  the  aid  of  the  above  notation,  nondimensional  variables  are  now  introduced  according  to 


(*i  ,  x2  ,  *3)  -  Cxi  ,  x2  ,  x3)/\  ,  t  -  t/tim)  , 


With  the  variables  defined  according  to  (5),  the  material  properties  are  seen  to  be  periodic  in  the 


*2,  x 3- plane  in  which  the  periodicity  of  the  fiber  lattice  structure  may  be  defined  by  the  cell.  It  is 


expected  that  stress  and  deformation  fields  will  vary  significantly  with  respect  to  two  basic  length  scales: 


(1)  a  "globaf  or  "macrtf  length  typical  of  the  body  size  or  loading  condition,  and  (2)  a  "micro"  length 


typical  of  "cell"  planar  dimensions.  Further,  it  is  expected  that  these  scales  will  differ  by  at  least  one 


order  of  magnitude  in  most  cases.  This  suggests  the  use  of  multivariable  asymptotic  techniques  (Ben- 


soussan,  Lion  and  Papanicolaou,  1978,  Hegemier,  Murakami  and  Maewal,  1979,  Sanchez-Palencia, 


1980).  This  approach  commences  by  introducing  new  independent  microvariables  according  to 


*;  =  v«  • 


Therefore,  all  field  variables  are  considered  to  be  functions  of  the  microvariables  x\  and  X3,  as  well  as 


the  macrovariables  x; ,  /  -  1-3: 


/Cx,,  Xj,  x3,  /)  -  /*Cx,,  x  j,  x3,  xj,  xj,  t  ; «) 


XUS 


-  7  - 


Spacial  derivatives  of  a  function  /  then  takes  the  form 


(Xk>x},t  ;c) 


(7b) 


where  9(  )/6jc,  =  0.  By  introducing  the  notation  (  ),.  =  8(  )/dx*  equation  (7b)  can  be  rewritten 
as: 


/,  “  /?,  +  1  fj-  ■  (7c) 

« 

In  the  sequel  /*  will  be  written  as  /  for  notational  simplicity. 

The  operations  (7),  when  applied  to  all  field  variables,  lead  to  the  following  "synthesized’  govern* 
ing  field  relations: 

(a)  Equations  of  motion 


<T 


JrJ 


P " 


(8) 


(b)  Constitutive  relations 

aft'  -  A<“>  Sy  eft'  +  *«">  ; 


(c)  Strain-displacement  relations 


,fa> 


fa) 


+  i#  +  -1-  (i#  +  «#)| 


(9) 


(10) 


(d)  Interface  continuity  conditions 

h<D_m.(2>  rjprjn  on^.  (H) 

At  this  point,  the  variation  of  field  variables  which  satisfy  the  periodicity  with  respect  tox*  is  assumed. 
According  to  this  condition  field  variables  take  equal  values  on  opposite  sides  of  the  cell  boundary.  The 
premise  allows  one  to  analyze  a  single  cell  in  an  effort  to  determine  the  distribution  of  any  field  variable 
with  respect  to  the  microcoordinates  x*.  The  x‘  -periodicity  condition  is  motivated  by  the  Floquet  and 
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Block  theorems  (Brillouin,  1946)  for  harmonic  wave  in  periodic  structures.  Certainly,  it  eliminates 
boundary  layer  effects.  However,  it  is  expected  to  provide  a  good  model  for  the  global  wave 
phenomena  in  fibrous  composites  with  periodic  microstructure. 

For  the  construction  of  a  mixture  model  it  is  convenient  to  cast  the  field  equations  in  a  variational 
form  by  using  the  Reissner  new  mixed  variational  principle  (Reissner,  1984).  In  the  Reissner  varia¬ 
tional  principle  the  variations  of  displacement,  strain  with  (10)  as  definition  and  transverse  stresses,  i.e., 
all  stress-components  except  aft’,  are  considered.  Thus,  it  is  convenient  to  rewrite  the  constitutive 
relation  (9)  in  terms  of  the  axial  strain  and  the  transverse  stresses: 
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+  +  j  «fe>.  -  2*}?>  ( •  •  •  ))J^;  dx\ 

+  /  (8u/2)  -  8i4(,))7r  +  87T(w,(2)  -  «,(l>)  ds’  dx,  dx2  dx3 

j?  ' 

“  /.f  J*  |  JE,  JJ  Sm/(o)(—  p^uffldx]  dxjJ  dx,  <ic2  <ix3 

+// ( ijs  8“-(a)  ft***;  ,  d3) 

where  4 fa)  denotes  the  x2 ,  x3 -domain  of  the  cell  occupied  by  material  a  (Fig.  2),  a  If*  is  used  for  the 

v 

approximate  transverse  stresses,  Tja)  denotes  the  traction  vector  on  the  surface  9  VT  where  the  traction 
is  specified,  ds  is  an  infinitesimal  line  element  on  «^,  and  dA  is  an  infinitesimal  surface  element  on 
the  boundary  of  V  :9  V.  In  (13)  basic  variables  are  the  displacement  «>>,  the  transverse  stresses  cr^] 

and  the  interface  traction  vector  Tt.  The  Euler-Lagrange  equations  of  (13)  include  (8a),  (11a),  (12), 
and 

77  =  oJf'Vj0  on  &  .  (14) 

The  above  variational  equation  (13)  furnishes  a  tool  with  which  a  mixture  model  can  be  obtained 
with  appropriate  trial  displacement  and  transverse  stress  fields.  The  basic  requirement  for  the  variables 
is  the  x* -periodicity  condition  on  the  cell  boundary  94  .  The  microstructural  variation  of  the  trial  func¬ 
tions  can  be  obtained  by  the  asymptotic  procedure  (Murakami,  Maewal  and  Hegemier,  1981). 


3.  Asymptotic  Analysis 

The  premise  that  the  composite  macrodimension  is  much  larger  than  the  microdimension, 
«  «  1>  and  the  form  of  scaled  equations  (8)  and  (10),  suggest  the  expansion  of  the  dependent  vari¬ 
ables  in  the  asymptotic  series: 

I H  1^**  (x* ,  X/ ,  I ,e )  ■»  L  e"  {u,(B) ,  (Ty (x^ ,  X/ ,  r ) 
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ir  (15)  is  substituted  into  (8)-(ll)  and  the  coefficients  of  different  powers  of  €  are  equated  to  zero,  a 
sequence  of  problems  defined  on  the  cell  is  obtained.  The  first  of  the  equations  in  this  sequence  fur¬ 
nishes 

,  aj%j.  -  0  .  (16) 

Equation  (16a)  implies  that  u,%]  is  independent  of  xj  and  yields  with  the  zero-th  order  expansion  of 
(11a): 

-  UtM  ( xk  ,  t)  .  (17) 

The  remaining  systems  of  equations  obtained  from  (8)-(10)  are,  for  «  >  0: 
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To  be  added  to  the  foregoing  are  the  interface  conditions  and  the  x'  -periodicity  conditions  for  n  ^  0: 

Htft-id?,  on  &  »  (21) 

«/<?)  and  <TjWn)vJ2)  are  x’- periodic  on  dA  (22) 

The  first  set  of  microboundary  value  problems  (MBVP’s)  for  crfy)  and  called  the  0(1) 
MBVP’s,  is  defined  by  (16b),  (18b),  (19)-(20),  (21b),  (22b)  with  n  -  0,  and  (21a),  (22a)  with  n  -  1. 
The  0(1)  MBVP’s  are  excited  by  Uno)J.  Similarly,  a  sequency  of  MBVP’s  is  defined  for  each  n  from 
(18M22).  With  appropriate  integrability  and  normalization  conditions,  higher  order  terms  may  be 
computed  by  solving  the  MBVP’s.  In  particular,  the  0(1)  MBVP’s  are  the  ones  solved  for  the  0(1) 
homogenization  theory  proposed  by  Bensoussan,  Lion  and  Papanicolaou  (1978)  and  Sanchez-Palencia 
(1980),  and,  also,  form  the  basis  of  the  mixture  theory  which  may  be  classified  as  an  0(c)  homogeniza¬ 
tion  theory.  The  asymptotic  approach  yields  the  microstructures  of  displacement  and  stress  fields  after 
solving  a  multitude  of  MBVP’s  which  are  complicated. 


'  In  order  to  use  the  approximate  solutions  of  the  MBVP’s  in  the  course  of  developing  a  mixture 

model,  and  to  ease  the  burden  of  solving  the  MBVP’s  exactly,  a  variational  procedure  was  adopted  by 
Murakami  (1985)  for  laminated  composites  with  the  help  of  the  Reissner  new  mixed  variational  princi- 
i  pie  (Reissner,  1984).  A  similar  approach  is  adopted  here  for  fibrous  composites.  To  obtain  the  lowest 

order  mixture  theory  by  using  (13),  it  is  necessary  to  obtain  trial  displacement  and  transverse  stresses 
to  0(c).  In  the  sequel,  the  trial  functions  are  obtained  for  a  hexagonal  cell  with  a  concentric  cylinders 
j  approximation  as  shown  in  Fig.  2.  In  Fig.  2,  (r,0)  are  micropolar  coordinates: 

r  “  +  •X'j2  ^  tan  0  =  x\/x\  ,  (23) 

by  which  r  =  1  constitutes  the  cell  boundary  and  r  =  denotes  the  interface  & .  The  quantities 

|  n(a)  indicate  the  volume  fraction  of  material  a  and  satisfy 

i 

n(\)  +  n(D  =  j  (24) 

In  terms  of  the  polar  coordinates  the  x '-periodicity  conditions  for  a  hexagonal  cell  with  the  concentric 

i 

cylinders  approximation  reduce  to  the  form: 

f(xrj9,t)**f(xk,r,Tr  +  9,()  at  r-1  .  (25) 


4.  Trial  Displacements  and  Transverse  Stresses 

The  0(1)  stress  and  0(e)  displacement  fields  are  obtained  by  solving  the  0(1)  MBVP’s  which  are 
defined  by  (16b),  (18b),  (19M22)  and  (24).  These  MBVP’s  are  excited  by  Uh0)j.  The  exact  solution 
of  w,faj  is  furnished  in  the  Appendix.  For  the  mixture  formulation  it  is  convenient  to  introduce  an 
0(e)  displacement  variable  which  represents  Uno)J  +  Uj(0),i  according  to: 

S,  (xk,t)  =  J  u,(a)UjU)ds'  -  -j  f  uft\)v}x)ds'  (26) 

where  A  (— tt)  is  the  area  of  the  cell.  Due  to  the  fact  that  w(ftj  is  excited  by  Ut(0)j  +  Uj(0)j  one 


obtains 


Equation  (27)  can  also  be  obtained  if  one  substitutes  the  exact  Mjfa)  in  the  Appendix  into  (26)  and 
eliminates  Uno)J.  To  render  the  analysis  tractable,  it  is  preferable  to  utilize  an  approximate  form  of 
the  exact  solution  for  M/If)) .  The  exact  solution  indicates  that  the  following  form  of  the  0(e)  displace¬ 
ment  yields  a  good  approximation: 

M,fa)  0*  i')  “  ■$  Of*  ,t)g(a)  ( r )  cos  9  +  $  Or*  ,t)g<o)  ( r )  sin  0  (28a) 

where 


g(,)(r)  =  ~-  ,  S(2)0>  =  -4-  (~r  +  1)  .  (28b) 

Anticipating  the  0(e2)  difference  of  the  average  of  M,<a)  on  A  (a\  equations  (17)  and  (27)  yield  the  fol¬ 
lowing  trial  displacement  field: 

w,(“)  Oc*  jcr ,/)  »  Ujia)  Or* ,/)  +  €  M/fa)  Gc*  jcj.  j)  (29) 

where  M/fa)  is  defined  by  (28).  Equations  (29)  and  (28)  indicate  that  the  mixture  displacement  vari¬ 
ables  are  U,(l\  t/,(2),  5,  and  5/  with  the  constraint  (27). 

By  using  (29)  in  (19)  with  n  -  0  and  considering  the  0(e2)  differences  of  the  average  transverse 
stresses,  the  0(1)  trial  stress  field  may  be  expressed  as: 


O’22(o) 

(a) 

t220c*,/) 

(a) 

R  _ 

cos  20 

cos  20 

sin  20 

^33(o) 
^  23(o) 

* 

T330f*,/) 
T230f *,*  ) 

Of*,/) 

cos  29 
0 

+  fJJ’Or  *,0 

-cos  20 
sin  29 

+  t$(xk,t) 

sin  20 
0 

(30a) 


O'SKo) 

(a) 

[^31  Of*  ,/) 

(a)  ( 

8rt  j  1  ... 

sin  20 

cos  29  1 

O’l2(o) 

[T 12  Of*  ,/  ) 

+  -jj-  (//?  o*  ,o 

cos  20 

+  tJPOiS) 

-sin  29  f 

In  order  to  define  the  0(e)  trial  stress  field  it  is  convenient  to  define  the  0(e)  stress  variable 
according  to 


p'u*j)s-hl  -i/  »*wu*‘ 


(31) 


If  one  integrates  (8a)  over  A(a)  and  utilizes  the  x’  -periodicity  condition,  one  obtains  the  mixture 
momentum  equations: 
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«fa)  <r^f  +  (-1>*+1  P,  -  nfa)  p<“) 
where  the  average  operation  is  defined  by 

Wad)/..  —  1  _  f  r  /-(a)/— 


(33) 


From  (32)  it  can  be  seen  that  P,  represents  an  interaction  body  force  between  the  two  constituents 
across  the  interface.  Also,  the  form  of  (32)  with  P,  defined  by  (31)  satisfies  the  integrability  condition 
adopted  by  the  0(1)  MBVP’s  for  cr^])  which  are  defined  by  (18H20)  with  appropriate  n’s  and  (31). 
As  an  0(e)  trial  stress  field  which  satisfies  (31b)  one  may  use  the  following  approximate  fields: 
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As  a  result,  the  trial  transverse  stresses  are  expressed  as: 

^  **/V)  +  e  O-ffll  )(**  y*lj) 

where  5-,%)  and  ary\\)  are  defined  by  (30)  and  (34),  respectively. 


5.  Mixture  Equations 

By  substituting  the  displacement  and  transverse  stress  trial  functions  defined  by  (29)  and  (3S), 
respectively,  into  the  Reissner  variational  equation  (13),  one  obtains  the  following  relations  as  the 
Euler-Lagrange  equations: 

(a)  Equations  of  motion 

+  (-1  )"+,P,  -  ,  /  —  1  —  3  ,  (36) 

Mj  j  +  ±  (<r^>  -  <r{}a)  +  R  £»)  -  /  S(,„  .  /  -  1,2  ,  (37a,b) 
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kd  +  jr  *9“'  ~  <T$a)  +  ~ 1  • 1  “  1*3 
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(37c,d) 


where 


t(Mu  ,  Mn)  =  -|  I  ff  cr^g^icose  ,  sine)  rfrj  <*r3*  , 
w  "  /4  <*-i  ^ 


/  S  i  A«pw  ,  *(,)  “  7  .  *(2)  -  -^T  (2+  w<2>  +  ^  ,n  "<U)  ; 


(b)  Constitutive  relations 
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(43) 


where  it  is  understood  that 
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The  remaining  constitutive  relations  associated  with  <r}f)  are  obtained  from  (12a);  the  results  are 
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The  associated  boundary  conditions  are  on  9  V 
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(50a) 


(50b) 


where 


7<«p)  =  JLjfjT  7^°)  dx\  dx\  , 

.(r#.fJ-4  1  frfc«fa*.sin9)*S*;  .  (51) 

A  o—  1 

Equations  (36)-(50)  and  the  initial  condition 

t/>>  ,  Utf  ,  $  ,  at  /  -  0  on  F  (52) 


define  a  well  posed  initial  boundary  value  problem  with  respect  time  t  and  the  macrocoordinates  xk . 


6.  Harmonic  Wave  Dispersion  Spectra 

In  an  attempt  to  test  the  accuracy  of  the  mixture  model,  the  phase  velocity  and  group  velocity 
spectra  of  the  mixture  theory  have  been  compared  with  available  experimental  data  for  time  harmonic 
waves.  For  the  comparison  harmonic  waves  which  are  propagating  at  an  arbitrary  angle  of  incidence  in 
a  full  space  of  the  following  form  are  considered: 

U/fl)  ,  U[2)  ,  UP  ,  Up  ,  UP  ,  UP  ,  \/ik  ,  Sjik  ,  2$//*  ,  SjikY 

—  exp  |/fc  (x !  cos  0  +  x2  sin  <t>  cos  0  +  x3  sin  <f>  sin  0)  —  /<u|  U'  (53) 

where 

If  -  (t/P  ,  t/P  ,  'uP  ,  UP  ,  U$"  ,  UP  ,  5,  ,  sz  ,  4  ,  5*  ,  *F  (54) 

•  j 

and  [  ]r  denotes  the  transpose  of  (  ].  In  (53)  U{a)  and  s,  are  constant  amplitudes,  k  denotes  the 
wave  number,  o>  represents  angular  frequency,  <t>  is  the  azimuth  measured  from  the  x{  axis,  and  0  is 
the  longitude;  the  direction  of  the  wave  propagation  may  be  best  represented  by  the  wave  vector  k: 

k  -  k[cos«l>  ,  sin  4  cos0  ,  sin  0  sin  0)r  .  (55) 


Substitution  of  (53)  into  (36)  and  (37),  which  are  expressed  by  the  displacement  variables  with 


(40) -(46),  yields  an  eigenvalue  problem  for  eu  of  the  form: 

IK1U-  (ca>)2lM]U  (56) 

t 

where  (AT  ]  and  [A/]  and  are  11  x  11  real  symmetric  matrices,  the  elements  of  which  are  functions  of 
the  mixture  constants  and  the  wave  vector.  Furthermore,  [A/]  is  a  diagonal  matrix.  Upon  calculation 
of  the  eigenvalue  f<o  for  a  given  ek,  one  obtains  the  phase  velocity  Q  as 


Cp  —  (ea>)/(ek)  . 

For  each  computed  eigen  pairs  (e<o  ,  U)k,  k  -  1,2,  •  •  •  11  the  group  velocity 

r  dto 
*  “  dk 

can  be  obtained  by  taking  the  derivative  of  (56)  with  respect  to  ek: 


(57) 


(58) 


\K\  [Tk  -  2(€<o)Q  [A 1)  +  (€«)2lA/1  £4 


(59) 


For  the  kth  eigenpair  equation  (59)  yields 


UT{[ *1-  («w)W)}£/ 

(Ct\ - ; - ; -  .  (60) 
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In  the  subsequent  simulation  a  typical  cell  dimension  A  was  chosen  to  be  a  cell  radius  by  introducing 
the  concentric  cylinders  approximation  of  the  equal  area.  The  reference  elastic  modulus  and  density 
used  for  the  scaling  are 

E(m)  “  i,  ,  P(m)  «  I  P ^  (6D 

O-I  0-1 

where  E(a)  is  Young’s  modulus.  The  dimensional  frequency  v(Ht)  can  be  computed  from  c<o  by 

V-  (<  W  K/  £(m  )/P  (m  >/  (2w A  )  .  (62) 


Numerical  results  are  presented  for  a  boron-epoxy  composite,  for  which  experimental  results  were 
presented  by  Tauchert  and  Guzelse  (1972)  for  a  waveguide  case  <£  -  0°  and  a  waveflect  case  <f>  —  90°. 
The  material  properties  are  summarized  in  Table  1  in  which  the  values  for  Poisson's  ratio  are 
estimated.  In  the  simulation  A  was  computed  from  the  fiber  diameter  (-  2V«(1)A)  which  was 
1.016  x  1CT4  m.  The  group  velocity  spectra  for  a  waveguide  case  <f>  -  0  «  0°  are  shown  in  Fig.  3  for 
two  acoustic  modes:  a  "gross?'  longitudinal  mode  and  a  "gros^  shear  mode.  In  the  figure  the  same  sym¬ 
bols  as  the  reference  of  Tauchert  and  Guzelse  are  used  for  the  experimental  data  points.  It  is  noted 
that  reasonable  agreement  is  achieved  for  the  waveguide  case  in  which  pronounced  dispersion  is 
observed.  The  group  velocity  spectra  for  a  wavereflect  case  <f>  -  90°,  6  -  0°  in  which  the  wave  vector  is 
normal  to  the  fiber  axis  are  shown  in  Fig.  4  with  the  experimental  data.  The  figure  includes  three 
acoustic  modes:  a  "grostf’  longitudinal  wave  (P-mode),  a  "grossf  vertically  polarized  shear  wave  (SV- 
mode),  and  a  "gross"  horizontally  polarized  shear  wave  (SH-mode).  The  sets  of  experimental  data 
correspond  to  the  "grossf  P-mode  and  the  "gros^  SH-mode.  It  is  noted  that  there  are  significant  devia¬ 
tions  from  the  "gros^  SH-mode,  but  the  overall  agreement  is  not  unsatisfactory  if  one  admits  the  scar¬ 
city  of  the  experimental  data  and  the  difficulties  associated  with  the  measurement  of  shear  wave  veloci¬ 
ties.  It  was  reported  by  Tauchert  and  Guzelse  (1972)  that  a  shear  wave  exhibited  extremely  high 
damping  of  the  pulse.  A  similar  observation  and  the  scatter  of  shear  wave  data  were  reported  by  Sachse 
(1974)  who  conducted  modulus  measurements  of  boron/epoxy  composites  by  using  pulse-echo  tech¬ 
niques.  He  concluded  that  "the  measurements  of  the  present  investigation  indicate  that  shear  waves 
propagating  along  and  across  fibers  in  the  composite  materials  tested  do  not  always  propagate  at  the 
same  speed." 


Sutherland  and  Lingle  (1972)  reported  phase  velocity  measurements  for  tungusten/aluminum 
composites  whose  material  properties  are  shown  in  Table  2.  The  equivalent  cell  radius  A  was  computed 
from  the  given  fiber  spacings  which  yield  the  area  of  a  typical  cell  A  (-  n  A2)  0.579  x  10" 6  m2.  Fig¬ 
ure  5  shows  the  phase  velocity  vs.  frequency  relation  for  the  "gros^  longitudinal  mode.  A  reasonable 
agreement  is  observed  between  the  experimental  data  and  the  theoretical  prediction. 
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7.  Effective  Moduli 

The  0(1)  homogenization  theory  which  yields  the  effective  moduli  of  composites  can  be  obtained 


by  taking  the  limit  of  e  — * 

0  and  equating  the  constituents'  displacements 

U,(l)  -  -  Uj 

(63) 

By  introducing  the  above  constraints,  equations  (36)  yield 

<rfJ-p{m)UKn 

(64) 

where 

<r<w)  -  I  ,  p(m)  -  I  . 

0—  1  a*  t 

(65) 

Equations  (37)  yield 

<r^>  -  <r2Vfl)  +  Rtf*  -  0  ,  i-l,2, 3 

(66a) 

a- j,2a)  —  <rtfa*  +  Rtf*  —  0  ,  /- 1,3  . 

(66b) 

By  eliminating  S,  by  (66),  equations  (65a),  (40)  and  (47)  with  (63)  furnish 

(m )  —  )  Jgtw ) 
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where 
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?(m)-  Wu  ,  U2J,  t/3,3,  f/2.3+  >  t/3.1  UU,UU+  U2AV 


(68) 


and  [£(m)]  is  the  effective  modulus  matrix  with  transverse  isotropy  due  to  the  concentric  cylinders 
approximation  and  is  defined  in  the  Appendix. 

The  formulas  for  the  effective  moduli  (B2)  are  assessed  by  comparing  the  results  with  the  experi¬ 
mental  data  reported  by  Datta  and  Ledbetter  (1983)  for  boron/aluminum  composites.  The  results  are 
shown  in  Table  3  in  which  the  moduli  computed  from  the  effective  stiffness  theories  for  the  square  cell 
by  Achenback  (1976)  and  for  the  hexagonal  cell  by  Hlavacek  (1976)  are  included  by  using  the  formulas 


«} 

& 


reported  by  Datta  and  Ledbetter  (1983).  The  comparison  has  revealed  that  ail  high-order  theories  yield 
almost  similar  results.  It  can  be  easily  shown  that  the  formulas  for  the  effective  moduli  yield  values 
which  fall  between  the  upper  and  the  lower  bounds  obtained  by  Hashin  and  Rosen  (1964)  for  fiber- 
reinforced  composites. 


8.  Concluding  Remarks 

An  asymptotic  mixture  theory  with  multiple  scales  was  applied  to  unidirectionally  fiber-reinforced 
elastic  composites  with  periodic  microstructure.  In  the  model  construction,  Reissner's  new  mixed  vari¬ 
ational  principle  was  applied  to  the  synthesized  fields  with  multivariable  field  representations.  In  order 
to  assess  the  accuracy  of  the  model  the  mixture  dispersion  spectra  were  compared  with  the  experimen¬ 
tal  data  obtained  for  the  boron/epoxy  composite  by  Tauchert  and  Guzelse  (1972)  and  for  the 
tungusten/aluminum  composite  by  Sutherland  and  Lingle  (1972). 

A  satisfactory  correlation  with  the  experimental  data  indicates  that  the  proposed  mixture  model 
furnishes  a  basic  tool  by  which  dynamic  responses  of  the  composite  structures  can  be  investigated. 
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Appendix 

A.  Exact  14ft)  of  the  0(1)  MBVP's: 

“  G  g(a>(r )  { (U Ko) ,2  4-  U 2(0,1 )  cos  0  +  W j(o),i  +  U i(o)j)sin  0)  , 
u$)  “  4  (^2<0),2  +  ^3(0),3  +  A  £/,«,).,  )*<“>  OOcOS  9 

+  flia)[{^fa>(r)cose  4  k^gJr'ir)  cos  30}  (£/2(0).2  -  Um,3) 

+  {«<o)  (r)sin  9  +  Kfa)  gfr'  (r)sin  30}  (U2{ o>,3  +  l^ou)] 

+  62(o)[{3(1  -  K<“)Wa)(r)cos0+  (1  +  K^gfa'ir) cos  30}  (Um,2- ~  Gmj) 

+  {3(1  -  Kw)ftw  (r)sin  0  4-0  +  k (a))g,,f})  (r)sin  30}  (UWh3  +  C/3(o)  2)]  , 

m3(?1)  ”  b„  (U 2(0)  2  +  U 3(0) ,3  +  A  i/i(0),i)^<o)(/')sin  0 

+  ^  [{*<“>  (r)sin  0  -  K^a^sin  30}  (-  Um,2  +  UM ou) 

+  {*fa)(r)cos  0  -  k^gfrKr) cos  30}  {Um*  +  Vmj] 

4-  {g^  (r) cos  0  -  Kfa)a}“)  (r) cos  30}  (U2(0), 3  4-  £/3(0),2 >] 

+  A2<o)|{3(1  -  »c<“,)a/<a>(r)sin  0  -  (1  +  Kfc,)^>(r) sin  30} (-  (/2(o),2+  t/swu) 

+  {3(1  -  Kfa))S/(a)  (f >COS  0  -  (1  +  Kfa))aH)  COS  30}  ( i/ 2(0, ,3  +  (/j(0,.2)  (Al) 

where 

4  =  {(A  +m>(2)  “  (A  +M)(1)}/(2 </,)  ,  rf,  S  X  (A+M^Vn^  +  m(2,/(«(1,«(2)) 

a"  1 

<7  =  —  (t*U)  —  fi(2))/ d 3  ,  </}=  I  n^/n^  +  ni2)/(niun{2))  , 

a—  1 

x  =  u(l)  -  a<2>)/{(a+m)(i)  -  (a+m)<2,1  . 


(\+^)<")/M<">  ,  KUI=  (1— K(,,,)/(l+K<")) 


(A2) 


£/U)(r)  -  gj>)Hr)  -  r3  ,  g/)U(r)  -  0 


gii2)(r)  -  -  r3  +  r~l  ,  g/j2)  (r)  -  (-  r  1  +  r~3)/n(2)  , 

rJ  —  (1+k ®)'2 {4(l-(c ®+k 0)2)/''} -  3(1— ic(2))2r-*}  ,  (A3) 

In  Eqs.  (Al)  a{a)  and  b^a)  are  obtained  by  solving  the  linear  equations  for  x  —  [a{x) ,  b{[) ,  a$2) , 

b\2)V: 

I  A]  x  *=  B  (A4) 

4x4" 

4x1  4x1 

where 

1  ,/I21-0,/I31-M<1>/fl<"  ,/l4.  =  0 
A\2mt  3(1-k(,>)«<,)J  ,  A21=  <1+k(»)/i(,)J 

A32--  A 42  -  3n  wn  <n(l— k  (1))  ,  A  i3—  1  ,  /la--Kfa)/nfa)  , 

><33-M(2,/»(J){1-x<“>/»<,))  ,  A 43—  3ki2)fii2)/ nil)1 

Aia-~  3/i<2,(1  +  w(1,)(1-k(2)) 

A 24  —  (1  +  k(2,)«<1,J  —  (4 -  3n(I))(l  —  k<2)  +  K(2)J)/{n(l)0  +  k<21)} 

>4  j4 **  3#i<2,(l  —  K(2,)(rt(l)—  K(2>/nn>) 

A44--  ^<2,(1-K(2»)(n("+  3k<2,/«<1,)+  lV(2,(l-K<2)+»c(2,J)/{n0)J(l  +  K(2)))  ,  (A5) 

and 

Bx-B2-Ba-  0  ,  B}  «-  VU)  -m(2)>/2  .  (A6) 


It  is  interesting  to  note  that  for  most  of  practical  composites  b\a\  a  -  1,2  are  small  compared  with 


B.  The  Definition  of  [E (m)]  in  (67) 


(m) 

E„ 

E\2 

En 

0 

0 

0 

(m) 

Uu 

o-  22 

En 

E22 

E22 

0 

0 

0 

^2.2 

*33 

En 

En 

E22 
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0 

t/3.3 

o-  23 
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0 

0 

0 

E  44 

0 

0 

i/2,3  +  i/3,2 

0 

0 

0 

0 

Ess 

0 

i/3.1+  t/l.3 

0 

0 

0 

0 

0 

Ess 

t/..2+  ^2,. 

where 

E{f  -  £  n<a)(k  +  W(a)  ~  (X(l)  -  X(2))2M  , 

a—  1 

Elf  -  I  /j<“>X<“)-  (X0)-X<2)){(X  +  m)<0-  (X  +  m)(2)J/^i  , 

a—1 

£i7»  -  £  /»fa>(x  +  l»ya)  -  Ux  +  m)(,)  -  (x  +  m),2)}2M  -  Wn  -  m'21)2/^  , 

a—  1 

£jf 2  *■  £  flfa)X <a)  —  {(x  +  /i)<n  —  (X  +  fi)<2>) 2/d\  +  0i.(I)  —  nl2))2/d2  , 

a- 1 

Elf  “  (£22°  ~  Elf)/ 2  , 

-  I  <^(,>-/i(2))2/rf3  ,  (B2) 


d2  -  X  +  (X+M)(2)/(2n(1>«<2)) 

a-  1 


(B3) 


and  where 
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Table  1.  Material  Properties  of  the  Boron/Epoxy  Composite  Tested  by  Tauchert  and  Guzelse  (1972) 


Volume 

Fraction  n<a) 

Young’s  Modulus 

£<“> 

Poisson's  Ratio 

„(a) 

Mass  Density 

-(a) 

P 

^  Boron 

0.54 

379.2  GPa 

(55  x  ltf  psi) 

0.18 

2682  kg/m3 

(251  x  10-6  lb  sec2///?4) 

^  Epoxy 

0.46 

5.033  GPa 

(0.73  x  1<?  psi) 

0.40 

1261  kg/m3 

(118  x  1(T6  lb  sec2///!4) 

Table  2.  Material  Properties  of  the  Tungsten/ Aluminum  Composite  Tested  by  Sutherland  and  Lingle 
(1972) 


Volume 

Young's  Modulus 

Poisson’s  Ratio 

Mass  Density 

Fraction  n (nl 

£<“> 

P{a) 

^Tungsten 

0.022 

398  GPa 

0.28 

19194  kg/m3 

(2) 

'  '  Aluminum 

0.978 

71.0  GPa 

0.34 

2700  kg/m3 

Table  3.  Comparison  of  Effective  Moduli  of  a  Boron/ Aluminum  Unidirectionally  Fiber-Reinforced 
Composite  in  Units  of  10u  N/m 2 


Data3  Mixture  Square  Cell3  Hexagonal  Cell3 

Model  Model  Model 


E{f  2.450  2.551  2.480  2.551 

E\f  1.825  1.868  1.856  1.872 

E{ If)  0.779  0.661  - —  0.661 

E(f  0.604  0.578  -  0.578 

Eft'  0.526  0.604  - —  0.606 

Elf  0.566  0.559  0.451  0.561 


3  After  Datta  and  Ledbetter  (1983) 
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1972) 
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Figure  4.  Group  velocity  spectra  of  vavoreflect.  r.'.cd 
for  a  boron/epoxy  composite (Tauchert  and 
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Figure  5.  A  phase  velocity  spectrum  of  a  longitudinal 
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ABSTRACT 

In  order  to  improve  the  accuracy  of  the  in-plane  response  of  the 
shear  deformable  laminated  composite  plate  theory,  a  new  laminated 
plate  theory  has  been  developed  based  upon  a  new  variational  principle 
proposed  by  Reissner  (1984).  The  improvement  is  achieved  by  including 
a  zigzag  shaped  C*  function  to  approximate  the  thickness  variation  of 
in-plane  displacements.  The  accuracy  of  this  theory  is  examined  by 
applying  it  to  a  problem  of  cylindrical  bending  of  laminated  plates  which 
has  been  solved  exactly  by  Pagano  (1970).  The  comparison  of  the  in¬ 
plane  response  with  the  exact  solutions  for  symmetric  3-ply  and  5-ply 
layers  has  demonstrated  that  the  new  theory  predicts  the  in-plane 
response  very  accurately  even  for  small  span-to-depth  ratios. 

INTRODUCTION 

The  advent  of  metal  matrix  composites  and  their  application  in  the 
form  of  laminated  plates  has  created  a  demand  for  the  development  of  a 
laminated  composite  plate  theory  in  which  each  layer  may  experience 
plastic  deformation  with  constraint  hardening,  (Dvorak  and  coworkers, 
1976,  1984),  and  transverse  cracking.  In  order  to  simulate  the  inelastic 
responses  of  each  layer,  it  can  be  assumed  that  in-plane  strains  are  pri¬ 
mary  quantities,  while  transverse  strains  are  secondary  quantities. 

In  his  series  of  papers  Pagano  (1970,  1972)  obtained  elasticity  solu¬ 
tions  for  bidirectional  composites  for  the  problem  of  cylindrical  bending 
and  for  the  problem  of  the  simply  supported  rectangular  plate.  Pagano 
showed  the  importance  of  the  transverse  shear  effect  for  the  prediction  of 
accurate  plate  deflections  and  the  necessity  of  improving  assumptions  for 
in-plane  displacements,  which  are  assumed  to  be  linear  across  the  thick¬ 
ness  of  the  plate  in  the  Kirchhoff  as  well  as  the  Reissner-Mindlin  plate 
theories.  Since  the  development  of  laminated  plate  theories  including 
the  effect  of  the  transverse  shear  by  Yang,  Norris  and  Slavs  Icy  (1966), 
and  Whitney  and  Pagano  (1970)  many  higher  order  laminated  plate 
theories  have  been  proposed.  Historical  accounts  of  such  efforts  may  be 
found  in  the  articles  by  Seide  (1980),  Bert  (1984)  and  Reddy  (1984). 
However,  only  a  few  attempts  have  been  made  to  improve  the  in-plane 
strain  response. 

It  has  been  stated  by  Whitney  (1972)  that  the  accuracy  of  classical 
laminated  plate-  theories  for  detetminint  in-plane  stresses  is  ngt  improved 
jty  the  inclusion  of  shear  deformation.  Among  continuum  plate  theories 
in  which  the  number  of  equilibrium  equations  does  not  increase  with  the 
number  of  layers,  an  improvement  of  the  in-plane  response  was 
attempted  by  including  the  cubic  variation  of  the  in-plane  displacements 
across  the  thickness  of  the  plate  (Hilderbrand,  Reissner  end  Thomas, 
1949,  and  Lo,  Christensen  and  Wu,  1977).  An  a  posterior  estimate  pro¬ 
posed  by  Whitney  (1972)  uses  the  transverse  shear  strain  obtained  by  the 


first-order  shear  deformation  theory,  (Whitney  and  Pagano,  1970),  to 
obtain  in-plane  displacements  by  integration.  This  estimate  yields  an 
accurate  prediction  of  the  in-plane  strains.  However,  the  drawback  of  the 
approach  is  that  the  stress  fields  do  not  satisfy  plate  equilibrium  equa¬ 
tions. 

In  order  to  facilitate  a  theory  which  accurately  predicts  in-plane 
response,  a  new  laminated  plate  theory  has  been  developed  with  the  help 
of  a  new  variational  principle,  (Reissner,  1984).  The  improvement  is 
achieved  by  including  a  zigzag-shaped  C*  function  to  approximate  the 
thickness  variation  of  the  in-plane  displacements.  The  advantage  of 
using  the  new  Reissner  variational  principle  is  that  it  automatically  yields 
the  transverse  shear  constitutive  relations  with  appropriate  shear  correc¬ 
tion  factors.  A  comparison  of  the  in-plane  displacements  and  stresses 
predicted  by  the  proposed  theory,  with  Pagano's  exact  solution  of 
laminated  plates  in  cylindrical  bending  indicates  that  inclusion  of  the  zig¬ 
zag  function  predicts  the  in-plane  responses  more  efficiently  than  the 
inclusion  of  smooth  nonlinear  functions. 

FORMULATION 

Consider  a  laminated  composite  plate  composed  of  N  orthotropic 
layers  whose  principal  axes  coincide  with  rectangular  Cartesian  coordi¬ 
nates  X|,  X]  and  Xj.  The  coordinate  system  is  selected  with  xj  normal  to 


Fig.  1  Coordinate  system  and  the  approximation  of  in-plane  displace¬ 
ments 


267 


a 


the  plane  defined  by  x3  end  x3.  For  noutional  convenience,  (  )tt>, 
A-  1,2,  •  •  •  Ji  will  denote  quantities  associated  with  the  Ac  -th  layer.  For 
a  plate  of  total  thickness  A,  the  thickness  of  the  A-th  layer  is  denoted  by 
»W,A  as  shown  in  Fig.  1,  in  which  the  volume  fraction  nul  satisfies 

X  -  1  .  (I) 

*-i 

Unless  otherwise  stated,  the  usual  Cartesian  indicia!  notation  will  be 
employed  where  Latin  indices  range  from  1  to  3  and  repeated  indices 
imply  the  summation  convention.  In  addition,  the  notation 
(  )_,  s  8(  )/dx,  will  be  employed. 

With  the  aid  of  the  foregoing  notation  the  governing  equations  for 
the  displacement  vector  u,u>  and  the  stress  tensor  ojf'  in  the  A-th  layer 
become: 

(a)  Equilibrium  equations 

0  .  a?'-*?'  (2) 

where  /,  is  body  force; 

(b)  Constitutive  equations  for  orthotropic  layers 

[e„  op* 

*”1  0  0  Cm  It,, 


Kti  Ca 


2eu  - 


where  CJ*1  and  Cj*1  are  elastic  moduli  in  which  and  are 

the  reduced  stiffness  of  Whitney  and  Pagano  (1970); 

(c)  Strain-Displacement  Relations 

*</**  -  <«/*'  +  Uyj>)/2  ;  (4) 

(d)  Interface  Continuity  Conditions: 

„_<*♦»  t  *  —  1,2... JV—  I  ;  (5) 

(e)  Upper  and  Lower  Surface  Stress  Conditions 

<ri''-7r  on  x}  —  A/2  (6a) 

TT  on  x,--  A/2  .  (6b) 

The  objective  of  the  subsequent  analysis  is  to  derive  a  laminated 
plate  theory  which  would  have  improved  in-plane  normal  strain 
approximation  across  the  thickness  and  also  the  effect  of  transverse  shear 
deformation. 

The  basis  of  a  new  laminated  plate  theory  is  facilitated  by  a  new 
mixed  variational  principle,  (Reissner,  1984),  applied  to  the  N -layered 
composite  plate  whose  middle  plane  occupies  a  domain  D  in  the  X|,  x3- 
plane. 


f/|-U( 


8e/f>  +  <r#>  +  a»e&>  <r,?> 


+  +  28eif>  t#>  +  28ejf>  rjf> 


-  6rjJ»  (lift*  -  e#>  (•■•))  +  +  u£>  -  2eJJ>  (■••)] 

+  (“ft*  +  “ft*  ~  2ejf'(  -)jJ  dt,jd*idxj 


-  gUfCxiJfj,- -j)Tr  dxidx2 


Xl-4 


where  9Dy  denotes  the  boundary  of  D  with  the  outward  normal  v, ,  on 

which  tractions  T,  are  specified,  and  is  the  xj-domain  occupied  by 
the  A-th  layer.  In  (7)  «*(— )  implies  the  appropriate  right  hand  side  of 
(3b).  For  an  arbitrary  variation  of  u,  and  r«  which  is  the  approximation 
of  Oi,,  which  satisfy  (S)  and  (6),  it  can  easily  be  shown  that  equations 
(2a)  and  (3b)  are  obtained  as  the  Euler-Lagrange  equations  of  (7).  In 
the  process,  equations  (3a)  and  (4)  are  considered  to  be  the  definitions 
of  w/f*.  oi}',  and  oft*. 

TRIAL  DISPLACEMENTS  AND  TRANSVERSE  STRESSES 

A  new  laminated  plate  theory  which  includes  the  effect  of 
transverse  shear  and  the  improved  in-plane  normal  strain  approximation 
is  facilitated  \,y  introducing  appropriate  trial  displacement  and  transverse 
stress  field*  <nn  the  Reissner  variational  equation  (7).  As  a  trial  dis¬ 
placement,  a  zigzag  in-plane  displacement  variation  across  the  thickness 
of  the  plate,  whose  amplitude  is  expressed  by  S/UiXj),  has  been 
included  in  addition  to  the  linear  variation  of  the  Reissner-Mindlin  plate 
theory,  as  depicted  in  Fig.  1: 

Hu,lxi^r2^j)  -  V,  Ui  JCj)  +  (x),xj)x, 


+  $&,.*,)<- l)*-^xj*>/- 1.2 


uj*'  (X|  Jtj^rj)  -  t/jCx^j) 


xj*>sx3-x#>  (9) 

is  a  local  xj-coordinale  system  with  its  origin  at  the  center  of  the  A-th 
layer,  xjg’ .  The  inclusion  of  the  zigzag  function  was  motivated  by  the 
displacement  microstrocture  of  laminated  composites  (Murakami, 
Maewal  and  Hegemier,  1981).  In  order  to  test  the  efTect  of  in-plane 
zigzag  displacement  without  additional  complications,  the  trial  functions 
(8)  are  chosen  to  yield  the  lowest  order  theory  As  a  possible  higher 
order  theory  one  may  add  the  zigzag  functions  to  the  displacement  trial 
functions  of  Lo,  Christensen  and  Wu  (1977). 

Transverse  shear  stresses  are  approximated  by 

Ti*)_  fltt,U|A})9,*)(*j)  +  r,<*-”(X|jrjW‘)(xj) 

+  7-/*)(x,.XjW‘,Ui)  ,  I-  1.2  (10a) 

rj}'-0  (10b) 
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In  (10)  r,1*-”  and  7)u>  are  the  values  of  r*  at  the  upper  and  lower  faces 
of  the  *  -th  layer,  respectively,  and 

/  t£’*3  •  «2> 

a<*> 

From  the  definitions  of  ffu>  and  (6)  one  obtains 

r/*>-  IT  .  T,w-  7T  •  «3) 

Due  to  the  approximation  for  uj*1  which  yields  e/f 1  -  0,  cr33  becomes  a 
reactive  stress.  As  a  result,  <r33  is  obtained  by  integrating  the  *j- 
equilibrium  equation. 

LAMINATED  PLATE  EQUATIONS 

By  substituting  (8)  and  (10)  into  (7)  with  (3a)  and  (4),  and 
integrating  by  parts,  one  obtains: 

(a)  Equilibrium  equations 

+  Nh,j  +  (7^  -  7T)  +  F?  -  0  .  /-  1,2,3  (14) 

M„.,  +  A/l,j-Arj,  +  h(7r+  7T)/2+/*-0,/-l,2  (15) 

Lu.i  +  Lv.i-  (V-  (-l)*7r)  +  //'-0,(-l,2  (16) 

where 

(Nu ,  My ,  Ly)  &  Z  /  (Ijc,)  . 

(-1)* -^-xj*')ai*'<4r3,/J-l,2  (17a) 

W*.  K»)  =  I,  /  «.  (-  D*  -STZh^dx,  ,  (17b) 

*-t  jfa,  n^'h 

(T.tf.tf 0-  J  /  |l,*j  , 

;  °7C) 

(b)  Constitutive  relations 

(?/*’  -  •—  (7;(*-'>  +  r,tt)) 


-  C,'*'  (f/jj  +  ♦,)  +  (-1)‘  |  sj  ,1-1,2,  *-l,2„..«8) 

T5*(~  4  7’‘‘l  +  ("37*7  + 


4  trl) 


--jL)a(*,/c<*)  +  G<*+"/e,*+,>) . 

I  -  1,2  ,  *  -  1,  •  •  •  JV-  1  (19) 

where,  in  (18)  and  (19),  no  sum  is  implied  on  /,  and 

(B(lc8»+#/ica>) ,  (20 

Q,{t)  and  r,’*’  can  be  obtained  in  terms  of  V3J  +  ♦,  and  5,  by 
(18),  (19)  and  (13).  As  a  result,  and  K a  are  expressed  by 
U3J  +  ♦,  and  5,  through  (17b);  such  expressions  automatically  include 


the  appropriate  shear  correction  factors  by  virtue  of  the  Reissner  varia¬ 
tional  principle. 

The  appropriate  boundary  conditions  to  the  equilibrium  equations 
(14-16)  are: 

U,  or  Nupt+ N**!  .  /- 1,2,3  Ola) 

♦,  or  +  ,#-1,2  (21b) 

$  or  turi+f»»j  ,  #—1,2  .  (21c) 

The  remaining  constitutive  relations  for  N„%  Mu  and  Ly,  ij”  10 
are  obtained  from  (17a),  (3a),  (4),  (8)  and  (10b);  the  results  are: 

11(*) 

£u£ii 
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N»~  (wu+^8..) 


+  *#»  (tu  +  9u)  , 


+  (-  n* 
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.(*)  L  |Sl,|l  ] 
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M31-  £  »“>*  Cjt>  jjr#>  (l#u+  Uu) 
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+  (-l)*-2^(Su+S,.1))  , 
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CYLINDRICAL  BENDING  OF  LAMINATED  PLATES 

In  order  to  examine  the  accuracy  of  the  new  laminated  plate  theory 
cylindrical  bending  of  bidirectional  plates  is  considered.  The  plate  is  sim¬ 
ply  supported  on  the  ends  *  -  0  and  /,  and  is  infinitely  long  in  the  xy 
direction.  The  prescribed  boundary  conditions  on  the  upper  and  lower 
faces  are  given  by 

7?-«.  sini^.TT-TT-TT-O.  US) 

The  boundary  conditions  for  the  simply  supported- ends  are  from  (21): 

N„- (/,-Af„-L„-0  atx,-0 ,/  .  (26) 

For  the  clarity  of  presentation  the  constitutive  relations  (22-24)  are 
cast  in  a  nondimensional  form,  in  which  h  and  reference  moduli  £r  are 
used: 


Nn/(hEr) 

An  fin  0 

Uu 

Mn/(h7ET) 

- 

Bn  Cm  C\ j 

Lti/lhET) 

0  C„  Cjj 

Su 

where 

AU  -  I,  «»*  (C,V/£r) .  B,}~  1  «.<*>(*#>/*><£, 'fVfr)  . 

C„-I  n<‘>( (*#'/* )l+  nu,*/12}(Cn>/£r)  . 


C„  -  I  (-  1)*  (»<‘,1/6)(£1,f,/£r) .  Ca  -  I  (r,^/3)(C^/£T)  (28) 


Equations  (18)  and  (19)  may  be  expressed  in  a  matrix  form,  respec¬ 
tively: 

Q  -  \A )  T  -  di  ((/,.,  +  ♦,)  +  djfSi/h)  (29) 

Ifiir-  IF]  Q 

(30) 

where 

Q-  10 1(" .  Q <?,w’F/(h£r) 

r-  (r,<,,,r/1,,-  -.rp-,,F/£r 

(31) 

and  [A  J,  UM  and  IF)  are  matrices  of  N  x  (N-l),  (N-I)  x  (N-l)  and 
(N-l)  xN,  respectively.  In  (31)  I  F  implies  the  transpose  of  I  ). 
Solving  (29)  and  (30)  for  Q  one  obtains 

r-  lbf-*  if)  q 

(32) 

Q- ✓((/,.,  +  *,)+  S,(Si/h) 

(33) 

where 

((/,.  /,)  -  (l/l  -  HSl~'|Fl)— 1  W„  dj) 

(34) 

and  (/]  is  an  N  x  N  identity  matrix.  Finally  substituting  (33)  into 
(17b),  one  obtains 

|NJ1/(h£r)j  [DUD,,!  It/,.,  +  tp|l 

1  KiJEt  j  “  [D2,DaJ  [sjh  J 

where 

Dn-fl  d,  ,  Dn-ej  d, : 

Dn  "  ■  4\  .  Du  “  fl  ■  4i 

and 

?,  -  II.  1, ...  IF 


ft-2h  1 /*«».  1/»(1>, ...  (-  l)w/nw,F 


(35) 


(36) 


(37) 


For  the  cylindrical  bending  problem  equations  (14H16)  reduce  to 


N„.,-0 

(38a) 

^u.i  +  «•  rin  •y  x,  —  0 

(38b) 

o 

1 

s? 

1 

* 

(38c) 

Lti.,-*Ji-0 

(38d) 

From  the  boundary  condition  (26)  for  N,,  one  concludes  that 

Nn  -  0  for  0  <  x  <  /  . 
By  using  (39)  one  can  rewrite  (27)  as  follows: 

(39) 

\Cn-Bh!An  C„| 

[**ul 

(40) 

[Lti/ (hET)  j  [  C|j  cj 

t  1 

The  constitutive  relations  (40),  (35)  and  the  form  of  forcing  terms  in 
(38)  suggest  the  following  form  of  the  displacements: 

UJh  -  fi)  sin  (w //)jt|  ,  -  4,  cos  (*/l)x, 


SJh  -  svpos  (w//)x|  (41) 

It  can  be  easily  shown  that  equations  (41)  satisfy  boundary  conditions 
(26).  The  equations  for  the  unknown  constants  fij,  4,  and  r,  are 
obtained  by  substituting  (41),  (40),  (35)  into  (38b-d).  The  result  is 


p1DnPDn  Pf>n 

pD n  Dn+pHCn—Bh/Au)  Dn+piCn 

♦i 

- 

0 

PD 2i  Du  +  pJCu  Djj+pJCjj 

0 

where 

pm*h/l  .  (42b) 


NUMERICAL  RESULTS 

The  accuracy  of  the  proposed  theory  is  assessed  by  considering  the 
cylindrical  bending  of  a  symmetric  three-ply  laminate  for  which  an  exact 
solution  has  been  obtained  by  Pagano  (1970).  The  material  properties 
are  for  0  degree  layers 

tff/ET  -  £,?/£r  -  25.0627 


CW/Et-CW/Et- 0.5  (43a) 

and  for  90  degree  layers 

C}?/Er  -  1.0025  .  C&'IEt  -  0.2  .  (43b) 


In  the  following  figures  the  results  are  presented  in  terms  of  the  nondi- 
mensional  quantities  adopted  by  Pagano  (1970): 


5| 


Et  «i(0-*i) 

*9. 
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lOOFr*’ 


(j.0) 


?I1- •7-»ii(4aj)  .  *»“  -rffiiW-Xj)  . 

»•  *  4« 

Sj  -  xjh  .  (44) 

In  the  various  curves,  the  solid  line  indicates  the  elasticity  solution 
(Exact),  while  the  results  by  the  proposed  theory  and  the  first  order 
shear  deformation  theory  (FSD)  are  shown  by  a  broken  line  and  a 
dashed-dotted  line,  respectively. 

For  a  symmetric  three-ply  laminate  (0/90/0)  with  layers  of  equal 
thickness,  the  plate  deflection  •  the  span-to-depth  ratio  relations  have 
been  plotted  in  Fig.  2.  The  reeults  of  the  first  order  shear  deformation 
theory  have  been  computed  by  setting  *|  -  0  in  (42).  Since  the  present 
theory  is  a  higher  order  theory  than  the  first  order  shear  deformation 
theory,  the  present  theory  predicts  E}  more  accurately  than  the  FSD. 
The  variation  of  the  in-plane  displacement  S(  acroes  the  thickness  of  the 


are  N  x  l  matrices. 


Fig.  3  In-pline  displacement  variation  through  the  thickness  for  Fig.  5  Transverse  shear  variation  through  the  thickness  for  l/h  -  4 

l/h  -  4  and  6 
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pute  has  been  compand  in  Fip.  3*  and  b  for  //A  -  4  and  -  «,  rmpec- 
tively.  The  results  show  that  the  poor  approximation  of  the  PSD  hat 
been  significantly  improved  by  the  praaent  theory  which  indudet  a 
zigzag-ihaped  function  for  the  variation  of  S|  acroas  the  thicknem.  Ihit 
accuracy  it  reflected  directly  in  the  variation  of  the  in-plane  normal  ttrtat 
'ii  e*_*hown  in  Figs  4a  and  b.  The  variation  of  the  tiantvene  thear 
ttrets  tru  for  l/h  -  4  hat  been  thown  in  Fig.  5.  in  which  the  curvet  for 
the  pretent  theory  and  the  exact  toiution  are  almost  identical 


5-ply 
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lammate  (0/90/0/90/0)  with  layer*  of  equal  thickness  hat  beta  con- 
aidered.  A  companion  of  the  plate  deflection  it  diown  in  F*.  6.  The 
oonetponding  results  for  the  in-plane  diapiacementt.  the  in-plane  normal 
atreaaea,  and  the  trantverae  ttreaaet  are  drown  in  Fig  7,  Figs  I  and 
Fig.  9,  respectively.  The  ratulta  reveal  that  the  present  theory  predicu 
the  piste  deformation  very  accurately  even  for  imall  l/h.  Furthermore 
it  can  be  teen  from  Fip.  3  wtd  7  that  the  variation  of  the  in-pine  dit- 
ptaoementa  acroaa  the  thickness  of  the  plate  ia  approximated  better  by  the 
Bgzag-type  function  in  (la)  than  by  or  (anna. 
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Fig.  6  Comparison  of  mid-surface  trantverae  displacements 
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Fip.  I  In-plane  normal  stress  variation  through  the  thickness  for 
//*  -  4  and  6 
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Fip.  7  In-piine  displacement  variation  through  the  thickness  for 
l/h  “  4  end  6 


Fig.  9  Trantverae  their  variation  through  the  thickness  for  l/h  -  4 


CONCLUSIONS 

A  laminated  eompoaite  plate  theory  which  accurately  predicts  in¬ 
plane  response  has  been  developed  besed  upon  the  Reissner  new  varia¬ 
tional  principle  (1984).  A  unique  aspect  of  the  theory  is  that  it  includes 
a  zigzag-shaped  C*  function  to  approximate  the  variation  of  the  in-plane 
displacement  over  the  thickness.  The  accuracy  of  the  theory  has  been 
examined  by  considering  the  cylindrical  bending  of  laminated  plates 
which  has  been  solved  exactly  by  Pagano  (1970).  The  comparison  of  the 
plate  deflection  and  the  in-plane  displacements  and  normal  streases  for 
symmetric  three-ply  and  live- ply  laminates  has  demonstrated  that  the 
new  theory  predicts  the  in-plane  response  very  accurately  even  for  small 
span-io-depth  ratios. 
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ABSTRACT 


In  order  to  improve  the  accuracy  of  the  in-plane  responses  of  the  shear  deformable  laminated 
composite  plate  theories,  a  new  high-order  laminated  plate  theory  was  developed  based  upon  Reissner's 
new  mixed  Variational  Principle  [91.  To  this  end,  a  zig-zag  shaped  C°  function  and  Legendre  polyno¬ 
mials  were  introduced  into  the  approximate  in-plane  displacement  distributions  across  the  plate  thick¬ 
ness.  The  accuracy  of  the  present  theory  was  examined  by  applying  it  to  the  cylindrical  bending  prob¬ 
lem  of  laminated  plates  which  had  been  solved  exactly  by  Pagano  [1].  A  comparison  with  the  exact 
solutions  obtained  for  several  symmetric  and  asymmetric  cross-ply  laminates  indicates  that  the  present 
theory  accurately  estimates  in-piane  responses,  even  for  small  span-to-thickness  ratios. 


1.  INTRODUCTION 


The  increasing  use  of  composite  materials  as  thick  laminates,  in  aerospace  engineering  and  in 
automotive  engineering,  has  clearly  demonstrated  the  need  for  the  development  of  new  theories  to 
efficiently  and  accurately  predict  the  behavior  of  such  structural  components.  The  intrinsic  hetero¬ 
geneity  and  anisotropy  of  these  composite  structures  as  evidenced  in  the  stacking  of  several  fibrous 
layers  and  in  the  high  discontinuity  in  material  properties  across  the  interfaces,  make  the  classical 
theories  of  plates  and  shells  inadequate. 

The  inspiration  and  guidelines  for  the  subsequent  attempts  have  stemmed  from  Pagano’s  works 
[1,2,3]  where  the  exact  elasticity  solutions  for  the  problems  of  cylindrical  bending  and  simply  supported 
rectangular  plates  were  given.  Pagano  showed  the  importance  of  incorporating  the  effect  of  transverse 
shear  deformations  in  order  to  accurately  estimate  the  plate  lateral  deflection  and  the  need  to  improve 
upon  the  thickness  variation  of  the  in-plane  displacements,  which  are  assumed  to  be  C1  linear  functions 
in  both  classical  plate  theory  (CPT)  and  Reissner-Mindlin  plate  theory  (FSD). 

The  first  attempt  to  develop  a  general  linear  laminated  plate  theory  is  credited  to  Yang,  Norris  and 
Stavsky  [4],  Their  theory  is  an  extension  of  the  Reissner-Mindlin  homogeneous  plate  theory  as  applied 
to  an  arbitrary  number  of  bonded  anisotropic  layers.  Whitney  and  Pagano  [5]  extended  Yang,  Norris 
and  Stavsky’s  work.  An  important  conclusion  drawn  from  their  analysis,  which  was  also  emphasized 
later  by  Whitney  [6],  is  that  the  inaccuracies  of  the  classical  plate  theory  at  low  span-to-thickness  ratios 
for  determining  in-plane  stresses  are  not  alleviated  by  the  introduction  of  shear  deformations.  Whitney 
[6]  obtained  in-plane  displacements  by  integrating  the  transverse  shear  strains  deduced  in  [5].  This 
resulted  in  a  higher  order  approximation  which  accurately  predicted  in-plane  strains,  but  the  resulting 
modified  stresses  did  not  necessarily  satisfy  the  original  plate  equilibrium  equations. 

Since  then,  other  high-order  laminated  plate  theories  have  been  proposed  that  account  for 
transverse  shear  strains.  Of  these,  the  Lo,  Christensen  and  Wu  [7]  and  the  Reddy  [8]  high-order 
models  have  served  as  the  foundation  for  the  present  theory.  In  their  paper  [7],  Lo,  Christensen  and 
Wu  used  appropriate  higher  order  terms  in  the  power  series  expansions  of  the  assumed  displacement 


field  which  was  proposed  by  Hildebrand,  Reissner  and  Thomas  [8].  On  the  other  hand,  Reddy  [9] 
imposed  the  condition  of  vanishing  transverse  shear  strains  on  the  top  and  bottom  surfaces  of  the  plate. 
However,  this  theory  does  not  satisfy  the  continuity  condition  of  transverse  shear  stresses  at  the  inter¬ 
faces. 

The  objective  of  the  present  paper  is  to  improve  the  approximation  of  in-plane  variables  in  lam¬ 
inated  plate  theories.  In-plane  displacements  and  bending  and  stretching  stresses  are  considered  pri¬ 
mary  quantities  in  any  approximate  laminated  plate  analysis;  transverse  stresses  are  only  of  secondary 
importance  since  they  are  an  order  of  magnitude  smaller  than  the  primary  bending  and  stretching 
stresses.  By  using  a  new  mixed  variational  principle  proposed  by  Reissner  [10],  the  present  theory  is  a 
high-order  model  which  improves  upon  existing  theories  by  including  in  the  assumed  in-plane  displace¬ 
ment  variations  across  the  plate  thickness:  1)  a  zig-zag  shaped  C°  function  as  detailed  by  Murakami 
[11];  and,  2)  Legendre  polynomials.  The  advantage  of  using  Reissner’s  new  mixed  variational  principle 
is  that  it  automatically  yields  the  appropriate  shear  correction  factors  for  the  transverse  shear  constitu¬ 
tive  equations.  Another  attractive  feature  of  the  proposed  theory  is  that  the  number  of  equations  to  be 
solved  is  not  increased  as  the  number  of  layers  becomes  larger  and  larger.  A  comparison  of  the  pro¬ 
posed  theory  with  Pagano’s  exact  elasticity  solution  for  symmetric  and  asymmetric  laminated  plates  in 
cylindrical  bending,  shows  that  in-plane  displacements  and  stresses  are  accurately  predicted  by  the  inclu¬ 
sion  of  the  zig-zag  shaped  function  and  the  Legendre  Polynomials. 
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2.  FORMULATION 


Consider  an  N -layer  laminated  composite  plate,  shown  in  Fig.  1,  with  principal  axes  coinciding 
with  a  Cartesian  coordinate  system  OclTx2r*3),  such  that  the  x3-axis  is  perpendicular  to  the  plane 
defined  by  xx  and  x2.  The  following  notation:  (  )ik\  k  -  1,2,.. ,N  will  designate  quantities  associated 
with  the  /c^-layer.  The  thickness  of  each  layer  is  n(k)h,  where  h  is  the  total  thickness  of  the  plate. 
The  volume  fractions  n(k)  satisfy  the  relation 

I  n<*>«  1  (1) 

*-i 


Unless  otherwise  specified,  the  usual  cartesian  indicial  notation  is  employed  where  latin  and  greek 
indices  range  from  1  to  3  and  1  to  2,  respectively.  Repeated  indices  imply  the  summation  convention 
and  (  )  ,  is  used  to  denote  partial  differentiation  with  respect  to  x, . 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement  vector  U/(k) 

th 

and  stress  tensor  o',}*’  associated  with  the  A:  -layer  are: 
a)  Equilibrium  Equations 

<rJ/j  +  fi(k)  -  0  ;  or,j*)  -  <rj)k)  (2) 


where  /,  are  the  body  forces; 


b)  Constitutive  Equations  For  Orthotropic  Layers 


<T|1 

(*) 

Cu 

Cl2 

0 

(*) 

e\\ 

<*) 

Cij/C^ 

O’ 22 

- 

c,  2 

C22 

0 

en 

+ 

C23/C33 

O’ 12 

0 

0 

2^12 

0 

*33 

<*> 

C13/C33  C23/C33 

0 

<*> 

*11 

(*) 

I/C33 

0 

0 

(*) 

0-33 

2*23 

“  - 

0  0 

0 

*22 

+ 

0 

1/C44 

0 

0-23 

2*3. 

0  0 

0 

2*12 

0 

0 

1/C55 

0’31 

(3a) 


(3b) 


where  Q  are  the  elastic  constants  and  Cy  (/  J  -  1,2,6)  represent  the  reduced  stiffnesses  introduced  by 
Whitney  and  Pagano  [5]; 
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c)  Strain-Displacement  Relations 


euk)  =  7  («.\/)  +  ujV]  ; 


d)  Interface  Continuity  Conditions 


^  o-]*) « o-Wr-n  ;  ; 


e)  Upper  and  Lower  Surface  Stress  Conditions 


<rly ’  -  T+  on  x3-  — 


<rtt°  -  Tr  on  x3--|  • 


The  objective  in  developing  a  new  laminated  plate  theory  is  twofold:  first,  to  improve  the  assumed 
variation  of  in-plane  displacements  through  the  thickness  of  the  plate  and  second,  to  include  the  effect 
of  transverse  shear  deformation.  In  order  to  carry  out  this  task,  Reissner’s  new  mixed  variational  prin¬ 
ciple  [10]  was  applied  to  the  N -layer  composite  plate  whose  middle  surface  occupies  a  domain  D  in  the 
xt,  x2-plane: 

JJ  jf  |sej*)<r(5*)  +  [u**J  +  u3(k>  -  2e£}  +  [u3(V  -  ej*1  (  ••)]8t#)J  dx3  dxx  dx2 

-  ff[l  f  Su,lk)  f,ik)  dx3\dxldx2+  f  I  f  t>Uj(k)  Tj  dx3  ds  (7) 

D  |  *  A(k)  |  SPj-  [  *  ^ «*> 

+  fj  |sm,0)  (xi^2,-|)7^  -  Su,(;v)(xltx2— -|)  T~  dxx  dx2 

V 

where  bDr  denotes  the  boundary  of  D  with  outward  normal  va  on  which  tractions  7]  are  specified  and 
A <k)  represents  the  x3-domain  occupied  by  the  A;1*1 -layer.  Also  e3i  (...)  implies  the  appropriate  right- 
hand  side  of  (3b).  Due  to  the  nature  of  Reissner's  mixed  variational  principle,  Eqs.  (3a)  are  taken  to 
be  the  definitions  of  o-ffi  used  in  connection  with  (7). 


3.  TRIAL  DISPLACEMENT  FIELD,  TRANSVERSE  AND  NORMAL  STRESSES 


The  high-order  laminated  plate  theory  which  takes  into  account  the  effect  of  transverse  shear 
strains,  is  obtained  by  including  the  Legendre  polynomials  of  order  n  -  1,2,3  with  respect  to  the  xy 
coordinate  to  a  zig-zag  in-plane  displacement  variation  of  amplitude  S,  (xlrx2)  across  the  plate  thickness. 

The  appropriate  trial  functions  used  in  connection  with  Reissner's  mixed  variational  principle 
Eq.  (7)  are  taken  to  be: 

a)  Trial  Displacement  Field 

Ui<k)(xlrX2S3)  “  UiOc,^  +  (4»I<*I*I>M>  +  S/(*lr*2)(-  D*  — (TjT  X3(,t) 

i  /r*'A 

(8) 

+  rX2WO+ 


2x 

where  £  =  —r^-  and  P„  (£)  are  the  Legendre  Polynomials  of  order  n .  It  is  also  understood  that  <ft%  =  0. 
h 

x3((c)  is  a  local  x3-coordinate  system  with  its  origin  at  the  center  x^1  of  the  k  th-layer,  i.e. 

x3<*>  =x3-x#>  .  (9) 


Eq.  (8)  may  be  regarded  as  a  superposition  of  a  zig-zag  function  and  the  cubic  variation  as  proposed  by 
Lo,  Christensen  and  Wu  [7],  with  the  exception  that  here  Legendre  polynomials  are  used  instead  of 
single  powers  in  x3; 

b)  Trial  Transverse  and  Normal  Stresses 


Ti*)(xlrxj,x3)  -  QP  0n»*2)F|(r)  +  &P  (xiS2)F2(z)  +  JP  (xlrx2'>F3(z) 

(10a) 

+  trj*  0(xi^C2)+  TP  Cx, )  +  [T’j*  n  (xi^c2)  —  T<k)  bc^2)]Fs(z)  ; 

t#>(X|,x2,x3)  “  Qik)  (x\*2)F\(z)  +  R$k)  (jc\jc2)F$(z)  +  Jjk)  (xiJc2)F3(z)  +  IP  (x{ jc2)F2(z) 

(10b) 

-i-  (r3(*-1)Cx,^2)+  TP  (xiS2)]F*(z)  +  [r3(*  1)Cxi^2)-  TP  ^^2))Ft(z) 


where 


F,(r) 


5 

nlk)h 


2lz*  — 


,  F2(z) 


-  30 
(nlk)h)2 


(4 r3-  z) 
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F}(z)  =  (k)f,\3  2 4  ^2+  4^  ’  ^4(2)  “  35r4  ^  *2  +  ^ 


Fj(z)  -  10z3  -  -j  z  ,  F6(z)  -*  ~('w(*)^ )2  “  14*3  +  *|  2) 


^■7(2)  -  1„ik^4  012z5_  40z3+  3z)  ,  Fg(z)-  126z5  —  35z3+  -y  z 


xj*'  1  1 

^  2  =  ^h 


Also  ,  ((?/**  ,  /?/*>  ,  7,ft) )  S  f  (1,  xj*>  ,  x3<*)2 )  T#>  dbfj 

>!<*> 


/i*>  =  /  xj*>3  T#>  (&3 


In  (10)  7"/*  n  and  7)tt)  are  the  values  of  r3l  at  the  top  and  bottom  surfaces  of  the  /r1*1  layer  respec¬ 


tively.  From  (6) 


7J(0)  -  Tt  and  pS)  -  T~  . 


The  degree  of  the  polynomials  F)(z),  /— 1—8,  appearing  in  (10a,b)  is  consistent  with  the  order  of  trun¬ 
cation  in  the  assumed  expansions  (8)  for  the  displacement  m/*\ 


4.  LAMINATED  PLATE  EQUATIONS 

Substituting  (8)  and  (10)  into  (7),  using  Gauss’  Theorem  and  the  orthogonality  property  of  the 
Legendre  polynomials  one  obtains: 

a)  Equilibrium  Equations: 

Tr  +  F?-Q  (14a) 


MalJ,  -  N3i  +  j  (7T  +  Tr)  +  F»-  0 


where 


1aia  -  K3i  -  (T+  -  (-l)v  Tr)  +  0 

Lai»  -  3 M3i  +  ^  (Tr  -Tr)+Ff--  0 

Pm~  (5Ii,  +  -j-V3o)+-^  (T:  +  Tr  )  +  F'«0 
Maf),  Zap ,  Lap,  Pap  ,v  » 

Fv  f-W  rZ  rl  pP  =  ^  J 

r>  '  r<  '  r<  '  r>  '  r>  *“>  ^*) 

fl  ^  a  fy\  /  .>1  2x3  /  h  \1  n  /» 


tl,  f  !*,<*),<-  D*  .  (f 


,  q)2  P2(c)  ,  (|)3/>3({)]  ^(1)  dx3  (15a,b) 


(N3i  ,  M3(  ,  K3i  ,  I3( )  =  l  1,  j  ?,({),(-  D*  ,  (j  )2  />,(<;)  rfx3  ;  (15c) 


b)  Constitutive  Equations: 

•  For  Transverse  Stresses 


8  Jlk)  n^h 

(n(k)h  )2  30 


+  (7'0(*-,)  +  r<k) )  -  j  A  /!<*>  e.ft)  [  U3jt  +  Va  +  Sa(-  D* 


+  h  rtglk)(^3a  +  3£a)  + 
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1  d  (k)  n(k)2h  i-p(k-\)  _  T'(Ar)  \  =  ^  „(*)3  r 

H  R°  40  Ua  1  120  n  Cl 


<*) 


*3*  +  3£a  +  S3„  (- 1)* 


+  3h  n}k)  (£ia  +  5 <f>a ) 


(16b) 


Qi 


(k)  _ 


-1^2 -  +  (f(k-l)  +  j'(fc)  \  _  _  }hL  „(k)3  Q (k)  (c  +  56  ) 

(n(/c)/7)2  12  y  a  a  '  40  “ 


(16c) 


1 

)  /o (* )  _ 

5  Jjk>  3  /?i*» 

1 

1 

/M*) 

'“'Q 

12  Qa 

3(rt<A:)/i  )2  1  1  n(k)h 

12 

»(*+!) 


5  /<*+,)  3 


3  (nft+u/»)2  1  n(k+"h 

h 


(16d) 


126 

For  Normal  Stresses 


-  *(k) 

n  +  8 


/=•(*) 


„(*)  +  ^tt+i) 


•  (*) 


■  (*+!) 


_(*+i) 

■p(k)  _  _2 _  y (Ac-*- 1) 


'(*+1) 


8Jj*>  ,  /.<*>/» 


(yjAc-0  +  73(*)  )  «  2b.  „(*)  C3(*) 


Oik) _ ZU. _ L 

(«<*>/»  )2  30  '  ‘J  '  5 


V3  +  S3(-l)k  -J^  +  3h  n0(k>  fj 


+  «<*> 


U  +  h  n0(k)  V  +  - y  (3n„(*)2  -  -J-  )|  +  -y 


<*» 


(17a) 


_L  » (*>  _  32  _  7\<*> )  =  _LL  h2  „ik)3  r<k)  t, 

A*3  5  n(k)2h3  +  140  ^3  J}  >  350  h  "  C33  *3 


+  /j2w<*)3 

1050 


^  +  (-  n* 


n(k)h 


(17b) 


Oi*>  -  -LLU-  + 

(n(*>A)2  12 


I  D  <*>  _ 

,  *3  2nik)2h3 


96 


n<,(k)t  +  yp  (5nJk)2  -±)4> 

-  n(*>3  [£+  5/tno(k,0] 

(17c) 

i*’>-'2688 

(17d) 

- 11  - 


■  11 

Qlk)  Qlk+" 

+  11 

R$k) 

/?3*+1> 

,_5_ 

y3(*)  y3(*+i) 

12 

c3<r  cr° 

+  2  h 

«(*>c3(i» 

W(*+1>C3<3*+1) 

3  h> 

n<*)2C3(|)  +  n(*+,)2C3(jt+,) 

-  70 

nk)  /j*+,) 

h 

(*) 

JL _  7\(fc— 1)  1  1 A 

„(*)  „<*+!> 

v*  rw 

T(it)  i  ^  T(ifc+1) 

ll/i3 

n(k)3Cti'  w(*+i)3  C3(|+1) 

“  18 

eg'  h  +  10 

cru 

+  Cg+|>  r3 

where  in  (16a, b,c)  and  (17a,b,c,d)  k  ranges  from  1  to  N  while  in  (16d)  and  (17e)  k  ranges  from  1  to 
CV-1) .  Also,  no  summation  on  a  is  implied  in  (16)  and 


Clk)  =  Sal  Cir  +  Sa2  cl v  ;  njk)  =  V/h 


£/ 

t/.,i 

^2.2 

¥ 

*2.2 

- 

■S’  i.i 

*^2,2 

£ 

fi.t 

€  2,2 

* 

0u 

02.2 

(18) 


(19) 


By  solving  (16)  and  (17),  Q,(kK  R,(k),  J^k\  /j**  and  7](*)  are  obtained  in  terms  of  £4,  5,  ,  and 

0„  and  their  derivatives.  As  a  result,  the  quantities  ,V3/,  Af3/,  if3/,  I3l  of  Eq.  (15c)  can  be  determined 
as  functions  of  these  displacement  variables.  Such  expressions  will  automatically  include  the  appropri¬ 
ate  shear  correction  factors  by  virtue  of  the  Reissner  mixed  variational  principle. 


The  equilibrium  equations  (14)  are  supplemented  with  the  following  suitable  boundary  conditions: 


specify 

V, 

or 

-V./  * 

(20a) 

specify 

*, 

or 

^ai  V a  * 

(20b) 

specify 

St 

or 

Zai  » 

(20c) 

specify 

f, 

or 

^*ai  » 

(20d) 

specify 

0a 

or 

Vfi  • 

(20e) 

The  remaining  constitutive  equations  for  ,V0/J,  MaP,  Za/},  LaP  and  Paf}  are  obtained  by  substituting 


(3a),  (4),  (8)  and  (10b)  into  (15a)  to  yield: 


where  ,V«  LVn  ,V22  iVl2]r,  U  —  [£/u  t/2i2C/ 1,2+  f/2il]r  with  analogous  expressions  for  M,  P, 

■  LVj .  are  3  x  3  matrices,  [C](*’  is  a  15  x  5  matrix  an  V1* ,  ...,  V*  are  1  x  4  vectors,  which 


are  given  in  the  Appendix. 
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5.  CYLINDRICAL  BENDING  OF  LAMINATED  PLATES 

In  order  to  test  the  accuracy  of  the  present  theory,  cylindrical  bending  of  composite  plates  under 
sinusoidal  loading  is  considered.  The  plate  is  simply  supported  at  the  ends  x(  —  0  and  /  and  is  infinitely 
long  in  the  .x-j-direction.  The  prescribed  boundary  conditions  on  the  top  and  bottom  surfaces  of  the 
plate  are: 

77  -  0  ,  Tf  -  q  sin  on  x3  -  (22a) 

rr«r3-*0  on  x3  =  -  -J  .  (22b) 

The  boundary  conditions  for  the  simply  supported  ends  are,  from  (20): 


at  X!  —  0  ,  / 

(23a) 

1 

1 

N 

1 

jr» 

1 

1 

o 

at  xx  -  0  ,  /  . 

(23b) 

Using  surface  boundary  conditions  (22),  the  equilibrium  equations  (14)  for  cylindrical  bending  reduce 
to: 


1 

o 

(24a) 

'Vi3,l  +  ?  sin  “  0 

(24b) 

o 

1 

> 

1 

1 

(24c) 

Ml3A  ‘V33  +  2  q  sin  / 1  -  0 

(24d) 

Zt  u  —  A31  —  0 

(24e) 

t3.i  —  K 33  -  q  sin  ——  —  0 

(24f) 

Z.11,1  -  3A/  3j  -  0 

(24g) 

.y-jA.I' 
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The  matrices  i4 1 ],  [C3]  and  vectors  X|,  k4  are  given  in  the  Appendix.  The  vector  equations 

(26a,b,c)  and  (27a,b,c,d)  have  N -components,  while  the  vector  equations  (26d)  and  (27e)  have 
( N -  1)  components.  Matrices  lA  j],  ....  IC3 ]  depend  on  the  volume  fractions  n{k)  and  elastic  con¬ 
stants  ’ .  CiV  and  Cj*’ ,  while  the  vectors  Xi, ...,  *4  contain  the  displacement  variables  Ux, ....  £3 . 

Eqs.  (26)  are  easily  solved  by  substituting  Qu  Rx  and  Jx  in  terms  of  Tx  from  (26a,b,c)  into 

(26d).  This  yields  a  new  equation  involving  Tx  only,  which  can  thus  be  solved  for  T\.  Then  by  back 

substitution  expressions  for  Qx,  Rx  and  -jjj  Jx  in  terms  of  kx,  X2  and  X3  are  obtained.  Proceeding  in 

a  similar  manner  with  (27a,b,c,d)  Q 3,  —  /?3,  73  and  73  in  terms  of  ki,  k2,  k3  and  k4  are  deter¬ 

mined.  These  expressions  are: 

2']  [<j[/]-H0,])  -(j[/]-2L4(?1])|  X,| 

"  (jl/1 -2U0,])  -  (|[/)-4U0,])  +\2lARi 

—  LBC?1  ]  (Ax  -  2X3)+  0/1-  l BRx])k2  (29b) 


]  X2  (29a) 


({[/]-  [AQi  1)  -  (-jl/1  -  2L4C?3])1  U  UR3) 

a  A  „  +  —  jtjD  1  (8k4—  5k2)  (30a) 

(|l/]-2L403l)  -  (-jC/1  -  4L4C?3))  M  3  l2L4/?3]J  - 


(ifi/i+fw/ij) 
'•IT  I'1  *!“■!' 


(-fy17^  [k2 

(‘ffl/1  +  ‘^I®i?3l)  l'4 


j  [io[5(?3]1 

“  3T  [l6[5(?3]| 


where  [/]  is  the  iV  x  N  identity  matrix  and 


[AQi)  l4*,]l  [14  ill  ,  | 

uai  ur/t.]  -  [0,1  Wl  17BM7W]  <-u 


with 


[TV j]  -  {a[tq3\[ax\  -  -|j  [7H3Hfi,]  +  (C3])-‘ 


(31c) 


[AQ,  ],  [5/?,  ]  are  N  x  N  matrices,  while  [TVt  ]  are  C V-  1)  x  CV-  1)  matrices.  By  inserting  (29) 

and  (30)  into  (15c)  and  (21)  the  appropriate  constitutive  relations  for  the  cylindrical  bending  problem 
in  terms  of  the  displacement  variables  U\, ...,  f3  and  their  derivatives  with  respect  toxi  are  obtained. 

The  form  of  the  dependence  on  the  displacement  variables  U j, ...,  £3  of  the  constitutive  equations 
thus  obtained  and  the  nature  of  the  applied  load  suggest  the  following  expressions  for  the  displace¬ 
ments: 
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Si 
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tx/h 

Si 

li/h 

d»,/A2 

(32) 


where  the """  quantities  are  nondimensional  by  definition.  It  is  easily  proven  that  the  boundary  condi¬ 
tions  (23)  are  satisfied  when  (32)  are  substituted  therein. 

Finally,  inserting  (32)  into  the  constitutive  equations  obtained  in  the  manner  described  above  and 
these  in  turn  into  the  equilibrium  equations  (24)  and  (25)  yields  a  system  of  nine  algebraic  equations 
with  the  nine  nondimensional  quantities  U\,  f3  as  unknowns.  This  system  is  conveniently  written 
in  matrix  form  as 


[5]  U~  F 


(33) 


where  U-  [Ux  Vi  5,  £|  U3  ¥3  S3  ^l7" 


(34a) 


f-  (0,$,0,-j?,  0,  -<7,0,-^?,0]r 


(34b) 


and  [B)  is  a  9  x  9  matrix. 


6.  NUMERICAL  RESULTS 


In  order  to  assess  the  accuracy  of  the  present  theory  the  problem  of  the  cylindrical  bending  of  an 
infinitely  long  strip  under  sinusoidal  loading  is  examined.  The  exact  elasticity  solution  has  been  given 
by  Pagano  [1],  where  a  three  layer  cross-ply  laminate  was  considered,  the  0°  layers  being  at  the  outer 
surfaces  of  the  laminate.  The  elastic  properties  are: 


for  the  0°  layers 


~  «  25.062657  , 
ET 


^13 

Et 


0.335570 


— -  -  1.071141 
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(35a) 


and  for  the  90°  layers 


~  -  1.002506 
Et 


-^■-0.271141 

Z-T 


~  -  1.071141 
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where  Er  is  a  reference  modulus. 

We  follow  Pagano's  [1]  nondimensionalization  and  write  the  displacements  and  stresses  in  the 

form 


u[k)  - 


Et_ 

<7 


(0^3  ) 


c3(fc) 


Ej_ 

Q 


u?'  ’  0) 


(36) 


-  <*) 
<7\\ 


^  o-/*’  (-J 


Also  *3  -  Y  '  S  ~  J 


(37) 


In  the  various  curves  the  solid  line  represents  the  exact  solution  while  the  results  of  the  present  theory 
are  shown  by  a  broken  line.  Also  shown,  for  comparison  purposes,  are  the  results  given  by  the  first 
order  zig-zag  model  (11)  and  Lo,  Christensen  and  Wu’s  high-order  theory  (LCW)  (7),  which  are 


represented  by  a  dashed-dotted  line  and  dotted  solid  line,  respectively.  Symmetric  3,  S  and  9-ply  lam¬ 
inates  and  asymmetric  4  and  8-ply  laminates  were  examined,  to  test  the  present  theory. 

For  a  symmetric  3-ply  laminate  (0/90/0)  with  layers  of  equal  thickness.  Table  1  shows  the  values 
of  the  central  deflection  «3  obtained  from  the  different  theories  for  a  span-to-thickness  ratio  5  of  4  and 
6.  As  observed  the  present  high-order  theory  correctly  predicts  the  central  deflection  w3  to  the  first  two 
decimal  digits,  while  the  first  order  zig-zag  model  gives  a  better  result  than  LCW.  The  variation  of  the 
in-plane  displacement  ux  across  the  plate  thickness  is  compared  in  Fig.  2a  for  5-4,  where  it  is  seen 
that  the  curves  for  the  present  theory  and  the  exact  solution  are  almost  identical.  This  improvement  is 
also  reflected  in  the  variation  of  the  in-plane  stress  across  the  plate  thickness,  as  shown  in  Fig.  2b. 
Very  close  agreement  is  found  between  Pagano’s  exact  solution  and  the  present  theory,  which  has 
improved  upon  Lo,  Christensen  and  Wu’s  high-order  theory,  especially  at  and  in  the  neighborhood  of 
the  interfaces. 

The  present  theory  was  next  tested  for  a  symmetric  5-ply  laminate  (0/90/0/90/0)  with  layers  of 
equal  thickness.  The  central  deflection  w3  for  span-to-thickness  ratio  5  of  4  and  6,  is  shown  in 
Table  1  where  close  agreement  with  the  exact  solution  is  observed.  The  variations  across  the  plate 
thickness  of  in-plane  variables  u{k)  and  <?/**  are  compared  in  Figs.  3  and  4.  The  curves  for  the  present 
high-order  theory  and  the  exact  solution  are  again  almost  identical.  In  particular,  it  is  seen  that  the 
present  theory  has  considerably  improved  upon  Lo,  Christensen  and  Wu’s  model  in  the  interior  layers 
of  the  plate. 

To  further  assess  the  accuracy  of  the  present  high-order  theory  the  more  difficult  case  of  a  sym¬ 
metric  9-layer  cross-ply  laminate  (0/90/0/90/0/90/0/90/0)  was  considered.  The  0°  layers  have  equal 
thickness  h/ 10  while  the  90°  layers  have  equal  thickness  h/ 8.  The  results  for  the  central  deflection  w3 
are  given  in  Table  1  for  5  —  4  and  6  where  again  close  agreement  with  the  exact  solution  is  observed. 
The  variations  across  the  plate  thickness  of  the  in-plane  displacement  ux  and  normal  stress  <ru  are 
shown  in  Figs.  5  and  6,  for  5  —  4  and  6  respectively.  There  the  discrepancies  between  the  first  order 
zig-zag  theory  and  the  exact  solution  are  more  pronounced  than  in  the  3-  and  5-layer  cases,  as 


expected.  However,  the  results  of  the  present  theory  are  still  very  good  when  compared  to  the  exact 
solution. 

Finally,  asymmetric  4  and  8  cross-ply  laminates,  with  layers  of  equal  thickness,  were  examined. 
The  present  theory  predicts  accurately  the  central  deflection  uy  These  results  are  given  in  Table  2  for 
span-to-thickness  ratio  S  of  4  and  6.  The  variation  across  the  plate  thickness  of  the  in-plane  displace¬ 
ment  u(k)  and  normal  stress  trff*  are  shown  in  Figs.  7,  8  and  9  for  S  —  4  and  6.  From  the  curves  for 
u{k),  it  is  seen  that  the  first-order  zig-zag  theory  deviates  significantly  from  the  exact  solution  at  the 
bottom  layer  of  the  plate.  On  the  other  hand,  the  discrepancies  between  LCW  and  the  exact  solution, 
for  both  u[k)  and  a[k>  are  more  pronounced  in  the  interior  layers  of  the  plate,  while  the  present  high- 
order  theory  is  still  in  very  good  agreement  with  the  exact  solution. 


7.  CONCLUSION 


A  high-order  laminated  plate  theory,  which  accurately  predicts  in-plane  responses  of  symmetric 
and  asymmetric  laminates,  was  developed  with  the  help  of  Reissner’s  new  mixed  variational  principle 
[10].  The  improvement  was  achieved  by  including  a  zig-zag  shaped  C°  function  in  the  in-plane  dis¬ 
placement  variations  across  the  plate  thickness,  as  proposed  by  Murakami  [11],  while  the  non-linear 
variation  is  accounted  for  by  using  Legendre  Polynomials.  The  accuracy  of  the  theory  was  examined 
for  the  case  of  cylindrical  bending  of  an  infinitely  long  strip  and  compared  with  the  exact  elasticity  solu¬ 
tion  given  by  Pagano  [1].  The  results  for  the  central  deflection  and  in-plane  displacements  and  normal 
stresses  for  several  symmetric  and  asymmetric  cross-ply  laminates  indicate  that  the  theory  very  accu¬ 
rately  predicts  these  in-plane  responses  even  for  small  span-to-thickness  ratios.  In  all  the  cases  con¬ 
sidered,  the  proposed  theory  gave  better  in-plane  responses  than  thr  Lo,  Christensen  and  Wu  high- 
order  theory,  especially  in  the  interior  layers  of  the  plate.  It  was  also  observed  that  for  symmetric  lam¬ 
inates,  the  first  order  zig-zag  model  [11]  predicts  more  accurately  the  central  deflection  than  the  Lo, 
Christensen  and  Wu  high-order  theory. 
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APPENDIX 


•  Matrices  [iVj, ....  IP^]  in  Eq.  (21): 
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The  (  )',  (  )"  and  (  )"'  quantities  can  be  obtained  from  (Al,2,3,4)  by  replacing  therein  C/f*  by  C#’ 
C$*  and  C$6*  respectively,  where  k  ranges  from  1  to  N. 

•  Matrix  [C](**  and  Vectors  Vs , ...,  V*3  in  Eq.  (21): 
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•  Matrices  \A\\,  ■■■■>  (Cj]  in  Eqs.  (26)  and  (27): 
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•  Vectors  Xb  ....  #c4  in  Eqs.  (26)  and  (27): 

X[  —  h  (C/3.i  +  'tri)  0\  +  $i  b\  +  A2('Pj<i+  3|i)C|  +  A3^3,i  d\  +  A30,  e, 

\2  —  /l2(^3  !+  3f  i)/i  +  hS 3 1  £i  +  /t3(f  3,1+  50i)ft 
X3  A3  (*3.1  +  50,)/, 

K,  -  /»£/,,,  a2  +  a3  +  Si  63  +  A2*3c2  +  A2¥u  c3  +  A3* , ,  d3  +  A40, ,  ft 
¥2*  +  A2'?,.,  /3  +  AS,,,  ft  +  A3?,.,  ft  +  A40,,,  ft 

?♦"  “  A  0i,i  s3 
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Os 


(*>  C#>  ,  b{k>  -  |(-  1)*  C#»  ,  c{k>  -  j  nJk)n(k*C&  ,  d[k)  -  j-(3n0(*)2-  |)n(*)C#> 

}>  ,  of  -  ■|n“lcjsl  .  4i“  -  -j(-  1)*  eft’  ,  c!*>  -  -jflfifCj}*  ,  cj‘>  -  -|»„,,,nl*lc;si 

Cl!1  ,  ef  -  j(5»,f J-  -|o»,‘>)n(‘>  C,'}1  ,/?'  -  «<*»  C#>  (A9) 

1“  -  C»>  .«!*>- .of  ««»C»> 


-26- 


Table  1  Central  Deflection  «3  for  Symmetric  Cross-Ply  Laminates  in  Cylindrical  Bending  Under 


FIGURE  8a: 


Thickness  variation  of  in-plane  displacement  ui * 
of^an  asymmetric  4-layer  cross-ply  laminate  for 


FIGURE  8b: 
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FIGURE  9a:  Thickness  variation  of  in-plane  displacement 

of  an  asymmetric  8-layer  cross-ply  laminate  for  S 
*»  4. 


FIGURE  9b:  Thickness  variation  of  normal  stress  of  an 

asymmetric  8-layer  cross-ply  laminate  for  S  »  4. 
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Abstract 

Based  upon  Reissner’s  new  mixed  variational  principle  (5)  for  displacements  and  some 
stresses,  a  new  high-order  laminated  plate  theory  was  developed  in  order  to  improve  the 
approximation  of  in-plane  variables  of  the  shear  deformable  laminated  composite  plate  theories. 
To  this  end,  a  zig-zag  shaped  C°  function  and  Legendre  polynomials  were  introduced  into  the 
assumed  in-plane  displacement  distributions  across  the  plate  thickness.  A  comparison  with  the 
exact  solution  obtained  for  symmetric  3-  and  5-layer  and  asymmetric  4-layer  cross-ply  laminates 
indicates  that  the  present  theory  provides  a  tool  with  which  one  may  efficiently  study  the 
extraordinary  skin  action  in  laminated  composite  plates. 


Introduction 

Sandwich  constructions  and  multilayered  thick  laminates  exhibit  under  combined  loads,  a 
very  different  mechanical  response  from  homogeneous  isotropic  plates.  One  such  effect  is  the 
stressed  skin  action  which  may  be  caused  by  either  the  drastic  change  in  material  properties  of 
the  laminae  or  by  the  slip  at  the  interface  of  adjacent  laminae.  In  order  to  facilitate  a  tool  with 
which  one  may  study  the  former  type  of  skin  action  laminated  composite  plate  theories  were 
closely  examined.  In  a  series  of  papers,  Pagano  (3,4)  derived  exact  elasticity  solutions  for 
bidirectional  composites  for  the  problems  of  cylindrical  bending  and  simply  supported  rectangu¬ 
lar  plates.  Pagano  showed  the  importance  of  incorporating  transverse  shear  deformation  effects 
for  the  accurate  estimation  of  plate  lateral  deflection.  Further,  he  showed  that  the  linear  varia¬ 
tion  of  in-plane  displacements  across  the  plate  thickness  adopted  in  the  classical  Kirchhoff 
(CPT)  and  Reissner-Mindlin  shear  deformable  (FSD)  plate  theories  may  not  be  appropriate  to 
simulate  stressed  skin  action  in  composite  laminates.  In  addition,  Whitney  and  Pagano  (7)  and 
Whitney  (8)  pointed  out  that  the  inaccuracies  of  CPT  at  low  span-to-thickness  ratios  for  deter¬ 
mining  in-plane  stresses  are  not  alleviated  by  the  introduction  of  shear  deformations. 

The  purpose  of  the  present  paper  is  to  point  out  the  differences  in  bending  responses  of 
laminated  composite  plates  when  several  displacement  approximations  are  introduced,  with 
emphasis  on  the  stressed  skin  action.  A  new  high-order  laminated  plate  theory,  with  a  first- 
order  theory  as  limiting  case  is  derived  using  Reissner’s  (5)  new  mixed  variational  principle.  It 
is  a  variational  principle  for  arbitrary  displacements  and  transverse  stresses  only,  in  which  the 
original  3-dimensional  stress-strain  relations  can  be  used.  The  unique  features  of  the  present 
theory  are  the  inclusion  in  the  assumed  in-plane  displacement  variations  across  the  plate  thick¬ 
ness  of:  l)a  zig-zag  shaped  C°-function;  and,  2)  Legendre  polynomials.  A  comparison  of  the 
present  theory  with  Pagano’s  (3)  exact  elasticity  solution  for  symmetric  3-  and  5-layer  and 
asymmetric  4-layer  cross-ply  laminates,  indicates  that  the  in-plane  responses  are  more  accu¬ 
rately  estimated  by  the  inclusion  of  the  zig-zag  shaped  function  and  the  Legendre  polynomials 
than  by  using  smooth,  C' -interpolation  functions  (CPT,  FSD,1). 
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Problem  Statement 

Consider  an  JV-layer  laminated  composite  plate  of  uniform  thickness  h,  as  shown  in 


Fig.  1.  Plate  Geomeetry  and  In-Plane  Trial  Displacement  Field 

Fig.  1.  A  Cartesian  coordinate  system  is  selected  such  that  the  middle  surface  of  the  plate 
occupies  a  domain  D  in  the  xu  x2-plane,  the  x3-axis  being  perpendicular  to  this  plane.  The 
notation  (  )(k),  k  -  1,2,  •  •  •  ,N  is  used  to  designate  quantities  associated  with  the  k  th-layer. 
The  thickness  of  each  layer  is  nik)h,  such  that  the  volume  fractions  n(k)  satisfy  the  relation 

I  1  .  (1) 

*-i 

Unless  otherwise  specified,  the  usual  Cartesian  indicial  notation  is  employed  where  latin  and 
greek  indices  range  from  1  to  3  and  1  to  2,  respectively.  Repeated  indices  imply  the  summa¬ 
tion  convention  and  (  ) ,  is  used  to  denote  partial  differentiation  with  respect  to  x, . 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement 
vector  Uj(k)  and  stress  tensor  <r}k)  associated  with  the  k  th-layer  are: 

a)  Equilibrium  Equations 

Vp]  +  f,(k>  -  0  ;  <r,Jk)  -  <t fjk)  (2) 

where  /,  are  the  body  forces; 

b)  Constitutive  Equations  For  Orthotropic  Layers 


O’!! 

<*) 

C\\ 

Cl2 

0 

(*) 

*11 

(*) 

C,yc3  3 

0-22 

- 

c  12 

^22 

0 

*22 

+ 

Cii/Cu 

<*12 

0 

0 

cw 

2e .2 

0 

(*) 


(3a) 


*33 

(*) 

c,yc33  c2yc33  o 

<*> 

*n 

(*) 

i/c33 

0 

0 

<*) 

°"33 

2®  23 

-  - 

0  0  0 

*22 

+ 

0 

1/C44 

0 

°"23 

2e3i 

0  0  0 

2e  12 

0 

0 

1/C55 

°’3I 

ik) 


(3b) 


where  C,j  are  the  elastic  constants  and  Ctj  (ij-  1,2,6)  represent  the  reduced  stiffnesses  intro¬ 
duced  by  Whitney  and  Pagano  (7); 
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c)  Strain-Displacement  Relations 


<<!*’  -  7  (*>?’  +  «,?’)  ; 

(4) 

d) 

Interface  Continuity  Conditions 

u,(*)-u/*+ »  ,  < r#)-<r#+1)  ;  k-l,2,..JV-l  ; 

(5) 

e) 

Upper  and  Lower  Surface  Stress  Conditions 

<4°  -  IT  on  x3  -  -| 

(6a) 

-  7T  on  x3  -  -  -j  . 

(6b) 

The  purpose  of  the  following  analysis  is  to  develop  a  laminated  plate  theory  that  will 
improve  the  assumed  variation  of  in-plane  displacements  through  the  plate  thickness  and  that 
will  also  account  for  transverse  shear  strains.  To  this  end,  Reissner's  mixed  variational  princi¬ 
ple  (S)  for  displacements  and  some  stresses  was  applied  to  the  N -layer  composite  plate: 

fj  1 1  ^  1 8eijk)<rl^k)  +  lu<kJ+u&)  -  (  •)]8r£,+  [«#>  (...)]Sr#> 


dx. 


dx{dx2[ 


i 


8u/*)  /,<*>  dx 3 


dxldx2  +  f 

I  f  8 «/*>  T,  dxy 

9D 

T  k  Au) 

ds 


(7) 


+ 


C*t,*2 -  8u,W) (xirX2-j) 


dx,  dx 2 


where  8DT  denotes  the  boundary  of  D  with  outward  normal  va  on  which  tractions  7)  are 
specified  and  Aik)  represents  the  xj-domain  occupied  by  the  k  th-layer.  Also  e3i  (...)  implies 
the  appropriate  right-hand  side  of  Eq.  3b.  Due  to  the  nature  of  Reissner’s  mixed  variational 
principle,  Eqs.  3a  are  taken  to  be  the  definitions  of  c rffi  used  in  connection  with  Eq.  7. 


Trial  Displacement  Field,  Transverse  and  Normal  Stresses 

The  present  laminated  plate  theory  which  accounts  for  transverse  shear  strains,  is  a  high- 
order  theory  obtained  by  superposing  a  zig-zag  in-plane  displacement  variation  of  amplitude 
S/(x1tx2)  across  the  plate  thickness  to  the  Legendre  polynomials  of  order  ft  —  1,  2,  3  in  the 
variable  xj  (see  Fig.  1).  A  first-order  theory  (2)  can  be  recovered  from  the  following  equations 
by  disregarding  the  underlined  terms  and  numbered  equations. 

The  suitable  trial  functions  to  be  used  in  Reissner’s  mixed  variational  principle  Eq.  7  are 
chosen  as: 


a)  Trial  Displacement  Field: 


<*> 

U\  ^1  $1  £l  01 

«20flr>f2^3) 

“ 

U2  ^2  S2  02 

U}(Xi^C2dCi) 

u>t±b.k  0 

1 


(-1)* 


n(k)h 


Xj 


.<*> 


(-|)2  p2(0 

(-j)3  Pj(() 


(8) 


j 
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2xi 

where  {  s  — —  and  P„  (()  are  the  Legendre  Polynomials  of  order  n.  It  is  also  understood  that 

ft 

Ui ,  S,,  and  0O  are  functions  of  xx  and  x2  only.  x/k>  is  a  local  ^-coordinate  system  with 

its  origin  at  the  center  of  the  k  th-layer,  i.e. 


x?'=x3-xW  ; 


(9) 
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0 

6 

0 
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j° 

35 

0 

-15/2 
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3/16 

lo 

0 

0 

3 

0 

-1/4 
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0 

10 

0 

-3/2 

0 

(0 

0 

0 

0 

1 

0 

36 

0 

-14 

0 

5/4 

0 

-112 

0 

40 

0 

-3 

0 

126 

0 

-35 

0 

15/8 

0 

_  x?' 
n^h 


1 

2 


z<2 


Also  ,  «?/*> ,  *,<*> ,  //*> )  s  f  (1,  Xj(k)  ,  x^)2 )  T dbt3 

/<*> 


/j*>  3  /  *j*’*  ,g> 


dx  3 


.<*> 


(10) 


b)  Trial  Transverse  and  Normal  Stresses: 

+  /]“  (xiJCitei&b)  +  [7,,tt“l>(x|,x2>  +  r/*)0c1^2)]F4(z) 

+  lTp-n(XiJcd  -  7;(k)0r,^2)H(8M  +  8i2)Fs(z)  +  8,3  Fg(z)l 

where  Hy  is  the  Kronecker  delta  and 

/»(k)A  F,(z) 

(/»(k,/02  F2(z)/30 
(n<k)A)3  F2(z)/105j 
F4(z) 

/j(z) 

(/i<*»A)2  F6(z)/105 
(n(k)A)4  F7(z)/315 
Fg(z) 


(11) 


(12a) 


(12b) 


In  Eq.  10  r/*-l)  and  7](k)  are  the  values  U  t3(  at  the  top  and  bottom  surfaces  of  the  k  th-layer 
respectively.  From  Eqs.  6 

T,l0)  -  7T  and  T,{N)  -  7T  •  (13) 


In  Eq.  11,  for  the  functions  Ft  (z),  /  —  1,  4,  5  two  rows  of  coefficients  appear:  the  upper 
one  corresponds  to  the  high-order  theory  while  the  lower  one  corresponds  to  the  first-order 
„  theory  which,  in  addition  to  disregarding  the  underlined  terms  in  Eq.  8,  is  obtained  by  setting 

r  if3 = o. 

For  each  of  the  theories  considered,  the  degree  of  the  polynomials  Ft(z),  i  —  1,  ...,  8 
appearing  in  Eq.  11  is  consistent  with  the  order  of  approximation  used  for  the  displacements 
U/(k>  in  Eq.  8. 


Laminated  Plate  Equations 

Substituting  Eqs.  8  and  10  into  Eq.  7,  using  Gauss’  Theorem  and  the  orthogonality  pro¬ 
perty  of  the  Legendre  polynomials  one  obtains: 

a)  Equilibrium  Equations 

Na^+Tf-TT  +  F?- 0  (14a) 
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where 


3(«/4 

-  *  w  +  3  (7T  +  7T)  +  -  0 

(14b) 

Zaiji  ' 

-  Ky  -  (r,+  -  (-1)*  7T)  +  -  0 

(14c) 

Lai* 

-  3A/j,  +  -J-  (7T  -  7T)  +  Ft  -  0 

(14d) 

Pfi*S 

-  (si*  +  -j  n*) + y  (n  +  n ) + K  -  0 

(14e) 

^<»|8  *  ^afi  *  ^a$  *  P*  ^  <• 

fi".  /?.  Fh 

11, 4  Pi({),(-  D*  .  (|)2  /»2({) ,  (-|)3Pj({)) 

Ml 
l/,“’|  ** 

(15a,b) 

UVj,  ,  M3l  ,  K3i  yZjj  tLy)  =  I,  /  1. 1  1)*  , 


(-»*  *ik)  ,  (-|  )2  P2(0  t#>  dx}  ; 


b)  Constitutive  Equations 
•  For  Transverse  Stresses 


+  «1  ^k)h  (ri*"0  +  T™  \-byh  »<*>  Cj"  [  ^  +  ^«  +  5a(-  D* 
+  h  njk\V 3„  +  3ga )  +  £  (3  /.,<*>*  -  \  )g3„  +  ^  (S  /.0(*)2  -  |  )0tt  j 

{ ***  ~  JlWL  (T«k~l)  -  T°k) > "  »tt)J  e-tt>  +  $3*  (-  l)k 


3**a>  (fj-+50a)] 


<?Jk)  -  ^  (7'a(*-1)  +  nk)  n{k)3  C<k)  ({}„  +  50a )  (16c) 


— I—  L  n(fc)  4.  .  S^‘»<r> _ Lg“f-  + _ 1 — [ajO<t+,) 

ei*>  .  3(i»(*)/o2  7h<*>/»  ej*+i)  l  2y" 

5  y(*+ 1)  3  £  (*+l) 

*T7n(°+l)/i)2  *  7  n“k+l)h 

- A  — yUt-o+fr  f  ff(k--  +  ■w<*>1>  1  Tik)+b 2  — *+l> 

£<*>  /a  3  £(*>  f<*+,)  /a  °2  £(*+l>  /o 


•  For  Normal  Stresses  (High-Order  Theory  Only) 


Qik)~~ffih)i  +  J!W  (7’j*-‘>  +  rj*>  )--^  V3+S3(-l)k  -^  +  3/rn„(*>*3 
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I k&UMhRXUT 


(17a) 


+  -£.«[*  +  *  ^  *  +  y  (3»0(t>2  -  J  >1  +  -y  [5  Hoik)3  ~  «att)  j* 

/  (if )  «  « 


I  *?’  -  jpm^T + ^ <r^"  -  ^  >  -  &  *’ c»' «» 

*  +  (-  D*  yy  5  +  3A  n™i  +  (5«0(*)2  -  ±  )* 


+  -LL  A2„<*>3 
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(?(*>  -  +  -yr  (7’^'1>  +  7’i‘))---^-  /»(‘>3  H  +  5A  »,<*>*  1 


T  *i‘>  - 


15  /,< 


(*) 


_  n™h 

2/i(*,2A3  T  96 


+  -2-^-  (rj*-0  -  rjfc)  )  - 


-  ll/i4 
2688 


* 


(17b) 

(17c) 

(17d) 


-  11  <?3tt)  .  Qlk+l) 

+  JI 

Rik) 

12  Cjf*  Cj|+1) 

+  2  h 

n<*+»Cft+n 

3A2 


yjk)  /,<*+!> 


-  70 


Uh3 


nk) 


1?+ » 


/*<*>2C#>  +  n<‘+»2cru 


«(*»Cif)  C#+1) 


_h_ 

18 


,(*) 


cr 

,«r+l) 

+  cF17 


rj*-l)  + 10 
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,(*> 


c#> 


«(‘+,)  rj*+l> 
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(17e) 


1  1 

where  in  Eqs.  16a,b,c  and  17a,b,c,d  k  ranges  from  1  to  N  while  in  Eqs.  16d  and  I7e,  k  ranges 
from  1  to  (iV—  1) .  Also,  no  summation  on  a  is  implied  in  Eq.  16  and 

ty  m  6aI  C#>  +  «a2  C#  >  ;  n}k)  =  x&'/h  .  (18) 

Also,  U  =  Ctf*  £/u  +  C#’  U j,2,  with  analogous  expressions  for  ¥,  S,  f  and  The 
coefficients  ax,  bx,  a2  b2  and  b3  appearing  in  Eq.  16  are  given  in  Table  1. 


Table  1  Coefficients  ax,  bx,  a2,  b2  and  b3  Appearing  in  Eqs.  16 


«i 

*1 

h 

b) 

High-Order  Theory 

1 

30 

2 

5 

1 

12 

1 

126 

4 

63 

First-Order  Theory 

1 

12 

a 

1 

10 

1 

30 

2 

15 

c)  Boundary  Conditions 


specify 

u, 

or 

Nal  *m  . 

(19a) 

specify 

V, 

or 

Mai  va  , 

(19b)i 

specify 

S, 

or 

Zal  va  , 

(19c)i 

specify 

(1 

or 

Lai  va  > 

(19d) 

specify 

or 

Pffa  Vff  . 

(19e) 
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It  should  be  pointed  out  that  for  the  first-order  theory  the  subscript  i  in  Eqs.  14b, c  and 
19b,c  will  assume  the  values  1  and  2  only. 

The  quantities  Nu,  Mv,  Kv,  Z},  and  I*,  obtained  by  inserting  Qp\  Rt(k\  J,(k\  /j*> 
and  T,ik)  from  Eqs.  16  and  17  into  Eqs.  10  and  ISc  will  automatically  include  the  appropriate 
shear  correction  factors  by  virtue  of  Reissner’s  mixed  variational  principle  (5).  The  remaining 
constitutive  equations  for  Naf},  Ma$,  Zap,  Lafi  and  Pafi  are  obtained  by  substituting  Eqs.  3a,  4,j 
8  and  10b  into  Eq.  ISa. 

Cylindrical  Bending  of  Laminated  Plates 

As  an  illustration  of  the  present  theory,  cylindrical  bending  under  sinusoidal  loading  of  an 
infinitely  long  strip  in  the  x2-direction  is  considered.  The  plate  is  simply  supported  at  the  ends 
xt  —  0  and  /.  The  prescribed  boundary  conditions  on  the  top  and  bottom  surfaces  of  the  plate 


irxi 

7?  -  0 , 7?  -  <7  sin 
77  -  77  “0 


on  *3  -  — 
on  Xi  -  -  -|  . 


The  boundary  conditions  for  the  simply  supported  ends  are,  from  Eqs.  19: 

U3  -  -  S3  -  £ 3  -  0  at  x,  -  0  ,  /  (21a) 

Nu- Mn- Zn- Ln- Pn-0  at  x,  —  0  ,  /  .  (21b) 

Numerical  Results 

In  order  to  assess  the  accuracy  of  the  present  theory,  the  cylindrical  bending  problem 
under  sinusoidal  loading  of  an  infinitely  long  strip  is  considered.  The  exact  elasticity  solution 
was  given  by  Pagano  (3)  where  a  3-layer  cross-ply  laminate  with  the  outer  layers  oriented  at  0° 
was  examined.  The  material  properties  are  for  the  0°  layers 

-  25.062657  ,  ^  -  0.335570  ,  -  1.071 141  ,  -  0.5  ;  (22a) 

Et  Et  Et  Et 

and  for  90°  layers 

Vil  _  1  nn-icni  ^13  _  n  int/11  ^33  _  ,  ml  1>(1  ^35  _  Q  2  /“M 


-z11  -  1.002506  ,  -rr1  -  0.271141  ,  -  1.071141 

Et  Et  et 


where  Et  is  a  reference  modulus.  Adopting  Pagano’s  (3)  nondimensionalization,  the  displace¬ 
ments  and  stresses  are  calculated  in  the  form: 


u{k)  — 


Ej_  uP®**)  -<*,  _  Et\  JOOAi  (±  Q) 

q  h  ’  U}  q  l*  *  K  2  ’ w 


-(*)  _  1  _(*)  (J.  r,  ■) 

<T|1  O-II  V  2  *  *3 ' 

Also  X3  *«  ——  ,  S  “  ~  (24) 

The  present  high-order  theory  and  the  corresponding  first-order  theory  are  compared  with 
Pagano’s  exact  solution  (3),  Lo,  Christensen  and  Wu’s  high-order  model  (LCW)  (1)  and 
Reissner-Mindlin  first-order  shear  deformable  theory  (FSD).  For  each  of  these  theories 
thickness  variations  of  in-plane  displacement  u[k)  and  normal  stress  <rik>  of  various  symmetric 
and  asymmetric  laminates  with  span-to-thickness  ratio  5  of  4  are  shown  in  Figs.  2,  3  and  4. 

For  symmetric  3-  and  5-layer  cross-ply  laminates,  with  layers  of  equal  thickness,  the 
results  of  the  present  high-order  theory  are  in  excellent  agreement  with  Pagano’s  exact  solu¬ 
tion,  as  can  be  observed  from  Figs.  2  and  3.  In  particular,  it  has  considerably  improved  upon 
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Lo,  Christensen  and  Wu’s  high-order  model  at  the  interfaces  and  in  the  interior  layers  of  the 
plate-  j.^so  the  Present  first-order  zig-zag  theory  gives  reasonably  good  results  even  for  the 
more  difficult  case  of  the  symmetric  S-layer  cross-ply  laminate 


Fig.  2a.  Thickness  Variation  of  In-Plane 
Displacement  for  N  —  3  and  5  —  4 


Fig.  2b.  Thickness  Variation  of  Normal 
Stresses  for  N  -  3  and  5  -  4 


Fig.  3a.  Thickness  Variation  of  In-Plane:  F|8-  3b.  Thickness  Variation  of  Normall 

Displacement  for  V  —  5  and  5—4  Stress  for  N  -  5  and  5  —  4 


Fig.  4a.  Thickness  Variation  of  In-Plane; 
displacement  for  N  -  4  and  5  -  4 
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For  an  asymmetric  4-layer  cross-ply  laminate,  with  layers  of  equal  thickness,  the  present 
first-order  theory  deviates  considerably  from  the  exact  solution  at  the  bottom  layer  of  the  plate, 
while  the  present  high-order  model  still  gives  very  good  results  as  shown  in  Fig.  4.  On  the 
other  hand,  the  discrepancies  between  LCW  and  the  exact  solution  for  both  u{k)  and  <?u  *  are 
more  pronounced  in  the  interior  layers  of  the  plate. 

In  all  the  cases  considered  FSD  gave  poor  results  when  compared  with  the  exact  solution. 

Finally,  the  improvement  of  the  present  theory  can  also  be  seen  in  the  results  for  the  cen¬ 
tral  deflection  u3  of  the  plate  as  shown  in  Table  2.  It  is  worth  noting  that  the  present  first- 
order,  zig-zag  theory  gives  closer  values  to  the  exact  solution  than  LCW  for  symmetric  lam¬ 
inates. 


Table  2  Central  Deflection  u3  of  Symmetric  and  Asymmetric  Cross-Ply  Laminates  for 


Number  of  Layers  N 

3 

4 

5 

Exact  Solution  (3) 

2.887 

4.181 

3.044 

Present  (High-Order) 

2.881 

4.105 

3.032 

Present  (First-Order) 

2.907 

3.316 

3.018 

LCW  (1) 

2.687 

3.587 

2.597 

FSD 

2.262 

3.088 

2.412 

Conclusion 

A  new  high-order  laminated  plate  theory,  which  accurately  predicts  in-plane  responses  of 
symmetric  and  asymmetric  laminates,  has  been  developed  using  Reissner’s  mixed  variational 
principle  (S).  The  improvement  was  accomplished  by  introducing  into  the  in-plane  displace¬ 
ment  variations  across  the  plate  thickness,  a  zig-zag  shaped  C°  function  and  the  Legendre  poly¬ 
nomials  of  order  n  -  1,  2,  3,  as  detailed  in  (6).  The  theory  was  tested  by  examining  the  prob¬ 
lem  of  cylindrical  bending  of  an  infinitely  long  strip,  whose  exact  solution  had  been  given  by 
Pagano  (3).  The  comparison  of  the  central  deflection  and  in-plane  displacements  and  normal 
stresses  for  symmetric  3-  and  S-layer  and  asymmetric  4-layer  cross-ply  laminates  with  available 
first-order  (FSD)  and  high-order  (LCW)  theories  has  shown  that  the  present  theory  very  accu¬ 
rately  predicts  the  in-plane  response  even  for  small  span-to-thickness  ratios.  Thus,  the  theory 
provides  a  tool  with  which  one  may  efficiently  study  the  extraordinary  skin  action  in  laminated 
composite  plates. 
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